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ABSTRACT 

Archival Chandra observations are used to study the X-ray emission associated with star formation 
in the central region of the nearby SAB(s)cd galaxy NGC 2403. The distribution of X-ray emission 
is compared to the morphology visible at other wavelengths using complementary Spitzer, GALEX, 
and ground-based Ha imagery. In general, the brightest extended X-ray emission is associated with 
Hii regions and to other star-forming structures but is more pervasive; existing also in regions devoid 
of strong Ha and UV emission. This X-ray emission has the spectral properties of diffuse hot gas 
{kT ~ 0.2 keV) whose likely origin is in gas shock-heated by stellar winds and supernovae with <20% 
coming from faint unresolved X-ray point sources. This hot gas may be slowly-cooling extra-planar 
remnants of past outflow events, or a disk component that either lingers after local star formation 
activity has ended or that has vented from active star-forming regions into a porous interstellar 
medium. 

Subject headings: galaxies: individual (NGC 2403) — galaxies: nuclei — galaxies: evolution — X-rays: 
galaxies 



1. INTRODUCTION 

Whether merger-induced collapse 

(iBarnes fc Hern quist" 1992: Cole et al."2000^. bar-driven 
inflow (iKormendy & Kennicutt ,2004 ) , or dynamical 
friction (jNoguc HIlooo) builds structure at the center of 
a particular galaxy, it is likely that star formation and 
(if present) central black hole growth will be strongly 
regulated by feedback from massive stars (and AGN 
activity) . 

Since much of the current central massive object 
growth and star formation is occurring in small, low- 
density, disk-dominated s pirals rather tha,n the massive 
but fuel-starved ellipticals ffleckman et al.ir2004[ ). nearby 
late-type galaxies are the ideal laboratories to view the 
growth of galactic structure in the current epoch. Here, 
we investigate star formation, feedback, and the growth 
of the central region of NGC 2403. We focus on the 
unique perspective enabled by X-rays in viewing the vi- 
olent and inherently high-energy phenomena associated 
with these dynamical processes and compare X-ray be- 
havior to that exhibited at other wavelengths. 

NGC 2403 is a SAB(s)cd galaxy in the M81 group of 
galax ies {D = 3.2 Mpc, l'=l kpc. iMadore fc FreedmanI 
|1991|) . It is the most massive, ^lO"'^" Mq, galaxy in 
the second largest of three M81 subgroups. There are 
7 dSph and di rr satellites of NG C 2403 known within 
this subgroup (IKarachentsev et a l. 2002) ranging in size 
from~3xl0^ to ~3xl0^ Mq (p^arachentsev et al.ll200l . 

As with many disk-dominat ed late-typ e spirals, 
NGC 2403 lacks a central bulge (|Kentl I1987D but does 
host a luminous c ompact nuclear star cluster first identi- 
fied in IR images (|Davidge fc Courteau|[200l . The mass 
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and luminosity of this cluster ar e comparable to t hose 
found in many late- type spirals (jBoker et al.ll200^ but 
it is old er and less compac t than typical for late-type 
galaxies (jYukita et al.ll2007l ). The dominan t age of stars 
withi n the nuclear star cluster is '^1.4 Gyr (jYukita et al.l 
|2007[ ). There are numerous younger asymptotic giant 
branch, red, and early-type supergi ants within the inner 
disk with stellar ages of 100 Myr (jPavidge fc Courteaul 
120021 ). Similarly, the Ha surface brightness near the 
center is lower than in surrounding regions while the 
youngest (^2-10 M yr) and most massive Hii regions 
(iDrissen et al.lll999D are 0.7 to 1.6 kpc from the cen- 
ter. iDavidge fc Courteaul ()2002n speculate that this age 
gradient may be due to the growth of a superbubble in 
the central region that has quenched star-formation in 
the central region while triggering activity further out 
through compression of surrounding gas during bubble 
expansion. In this picture, nuclear star formation may 
be an episodic phenomenon with roughly a few 100 Myr 
interval between major star-forming events. 

Such a star form ation cycle require s a replenishing 
source of cold gas. iSheth et all ([2005h . using molecu- 
lar gas maps, showed that bulges can be built by bar- 
driven gas inflow but that the process requires of or- 
der a Hubble time in galaxies like NGC 2403. Other 
studies have shown that the molecular and atomic gas in 
the central regions of NGC 2403 am ount to only a small 
fracti on of the dynamical mass (e.g.. lThornlev fc WilsonI 
Il995f ). I mportantly, the central total gas surfac e den- 
sity ( e.g., lThornlev fc Wilson 1995: Martin fc KennicuttI 
120011 ) ■ at least on kiloparsec sp atial scales, is be low the 
critical value for star formation (jKennicuttlll989l ) yet the 
global star-formation rate in NGC 2403 is a moderate 
1.2 Mq yr-\. 

More pertinent to the current study are several in- 
dependent investiga tions of the d ynamics of the Hi 
gas in NGC 2403 (ISickintd [l997t iSchaap et all [2000l : 
iFraternali et al.l l2002b[ ). The Hi rotation curve shows 
a flat gravitational potential typical of late-type spirals 
but there is also slower-rotating neutral hydrogen ex- 



Fig. 1.— Left: The GALEX FUV image of NGC 2403. The 21.'9 Xl2f3 eUipes makes D25 isophote of NGC 2403. The central 6'x6' 
region is shown as a square. The gray scale indicates intensity level in the unit of ct pixel"^. One pixel corresponds to I'.'S. Right: 
The same region in the merged X-ray data. The outer boundaries of the ACIS S2-S3-S4 combination are shown as thin-lined rectangles 
for each of the 4 overlaid observations. The central 6'x6' region is shown as a thick-lined square. The large circle in the center is the 2'.5 
radius circle used for the X-ray residual emission spectral analysis (see text). The two smaller circles and the small rectangle identify the 
background regions used for analysis of this residual emission. The gray scale indicates intensity levels in units of ct pixel"^. One pixel 
corresponds to 0'.'492. 



tending up to 15 kpc above the disk. Many properties 
of this an omalous Hi gas can be explained as galactic 
fount ains (jShapiro fc Fieldlll976t rBregmanlll980l : iSpitzeH 
Il990f ) but an additional, external, source of cold gas 
such as inf ailing clouds or even small satellite galax- 
ies might also be required (Fraternali fc BinnevI 120061 . 
[1)08; Struck & Smith 2009). High-resolut ion VL A ob- 
servations of NGC 2403 (|Fraternah et all |2002b) show 
very anomalous kinematic features on small scales that 
may be evidence of such clouds. Direct evidence of 
fountain activity is less conclusive. Measurements of 
Ha emi ssion line widths at se veral locations across the 
galaxy (jPraternali et al.l 12004 !) found only a few small- 
scale fountain-like features. To match the observed Hi in- 
flow with hot gas outflow would require many thousands 
of these small-scale fountains. Analysis of the spatial dis- 
tribution of diffuse hot X-ray-emitting gas in NGC 2403 
([Fraternali et al.ll2002af l concluded that as little as a few 
percent or as much as the majority of the inflowing Hi can 
be matched by outflowing diffuse hot gas, depending on 
assumptions made. Thus, it is not yet clear whether the 
buildup of the central regions of NGC 2403 is being fu- 
eled cyclically through galactic fountains, or by accreting 
gas from the intergalactic medium, or through mergers 
of small-scale gaseous satellites. 

If galactic fountains are at work and are the source of 
the majority of the observed infalling Hi gas, then there 
should be an imprint of this process in the X-ray emission 
from the disk correlated with the source of the fountains; 
namely localized star-forming regions that heat the gas 
through massive stars to the point of breakout from the 
disk. 

We confine our study of NGC 2403 to the inner 6' x 6' 
nuclear region (corresponding to 5.6x5.6 kpc; Figure [1] 
left panel) . This region is small enough to be fully imaged 
within a single observation by most instruments used in 
this study yet includes a large portion of the actively 



star-forming regions of NGC 2403 including 5 of the 6 
giant Hii regions identified bv iDrissen et al.l ([l999). 

We begin (§ ^ with an independent analysis of the X- 
ray data; examining X-ray-detected discrete sources and 
the underlying residual X-ra ^ emission. This analysis 
extends the previous work of Schlegel fc Pannuti[ ([200I 
and of Fraternali et al. (2002a) by including 3 subsequent 
Chandra observations which allow us to better quantify 
the spectral and temporal behavior of the X-ray emis- 
sion. Our results are consistent with both these previ- 
ous investigations. We expand our analysis to include 
optically-identified supernova remnants (SNRs) and Hii 
regions. Although individually X-ray-faint, we are able 
to characterize their bulk (average) X-ray temperatures 
and other X-ray properties using stacking analysis. 

We then turn to analysis of individual massive star- 
forming regions (§ [3|) in an attempt to determine the 
X-ray emission properties of these regions as a function 
of their age, mass, and extinction properties. We use 
a combination of (ground-based) Ha, mid-IR (Spitzer), 
and UV (GALEX) images to define these regions and to 
determine their basic physical characteristics. We then 
use our knowledge of the differences in X-ray properties 
between Hii regions and SNRs assembled in § [2| to inter- 
pret the properties of the hot gas. 

2. Chandra X-RAY OBSERVATIONS AND PRELIMINARY 
ANALYSIS 

In this section, archival Chandra images are exam- 
ined to derive the X-ray properties of four different 
source populations in NGC 2403; namely, bright X-ray- 
detected sources (mainly luminous X-ray binaries and 
background AGNs), optically-identified SNRs, optically- 
identified Hii regions, and the underlying unresolved X- 
ray emission. This allows us to parameterize these dif- 
ferent types of sources in terms of their X-ray temper- 
atures, luminosity distributions, and emission measures. 
This information will be used in § [3| to investigate in- 
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TABLE 1 
Chandra Observation log 



Date 


ObsID 


Instruments & Mode 


CTI 


2001/04/17 


2014 


ACIS-S TE 


36.1 


ks 


2004/08/13 


4628 


ACIS-S TE 


47.1 


ks 


2004/10/03 


4629 


ACIS-S TE 


45.1 


ks 


2004/12/22 


4630 


ACIS-S TE 


50.6 


ks 



dividual massive star-forming regions defined at other 
wavelengths. 

NGC 2403 was observed in full-frame mode with Chan- 
dra ACIS-S on four occasions for a total of ^^180 ks 
(Table [T]). We obtained level 1 event lists for all four 
observations from the Chandra data archive*^ and re- 
processed them using the CIAO (version 3.3.0.1) tool 
acis_process_events and calibration database CALDB 
3.2.1. Reprocessing removed pixel randomization, ap- 
plied CTI- and time-dependent gain corrections, and re- 
moved events with bad grades or bad status bits as well 
as bad and hot pixels. We created lightcurves for each 
observation using a 1 ks binning to check for periods of 
high background. We excluded intervals with total count 
rates >3cr above the mean rate for each observation. The 
final Good Time Intervals for the observations are listed 
in Table m 

SN 2004dj was used to define a common registration 
among the data sets. For this purpose, a circular Gaus- 
sian was fit to the image of the supernova to obtain 
an accurate centroid then the coordinates were adjusted 
to agree with the known position of SN 20 04dj (J2000 
^^^^jh^-jm^j Q^^s^ Decl.=-H65°35'57.84", lArgo et all 
I2004D . To register pre-supernova data, we bootstrapped 
using a bright point source near the aimpoint that is com- 
mon to both pre- and post-supernova data. The four 
co-aligned Chandra data sets were combined using the 
FTOOL utility fmerge to form a fifth merged dataset. 
The right panel of Figure [1] displays the merged X-ray 
image. The longest cumulative exposure is near the cen- 
ter of the galaxy where we have selected a 6'x6' region 
for analysis as indicated. Also shown are the galaxy's 
D25 ellipse, corresponding to the 25 mag-sec~^ contour 
in B. The 6'x6' region contains many point-like sources, 
extended star forming regions, and an extended region 
of X-ray emission near the galactic center. With the 
exception of the bright point-like X-ray sources, we are 
primarily concerned with the soft X-ray emission compo- 
nents of NGC 2403. For the soft emission, we fit spectra 
in the 0.4—2.0 keV range but quote results such as lu- 
minosities in the more familiar 0.5—2.0 keV band. For 
point sources, the fitting range is extended to higher X- 
ray energies. 



TABLE 2 
X-RAY Point sources 



R.A. 


Decl 


Count Rate " 




Variability 


J2000 


J2000 


10-* cts s 


lO^'^erg s-i 




07 36 24.5 


-1-65 


37 


13.1 


1.2 ± 0.7 


0.5 ± 0.3 




07 36 24.5 


-1-65 


37 


34.9 


-0.3 ± 0.4 


-0.1 ± 0.1 




07 36 25.2 


-1-65 


38 


40.6 


7.9 ± 1.7 


3.1 ± 0.7 


S 


07 36 25.6 


-1-65 


35 


39.7 


1158 ± 12 


456.6 ± 4.7 


L S 


07 36 25.9 


-1-65 


37 


57.0 


0.1 ± 0.6 


0.0 ± 0.2 




07 36 26.7 


+65 


36 


20.9 


0.8 ± 0.4 


0.3 ± 0.2 




07 36 32.6 


+65 


34 


52.6 


0.9 ± 0.3 


0.4 ± 0.1 




07 36 32.5 


+65 


38 


20.0 


0.1 ± 0.6 


0.0 ± 0.2 




07 36 32.5 


+65 


39 


00.7 


29.7 ± 3.0 


11.7 ± 1.2 


S 


07 36 33.0 


+65 


34 


57.8 


3.9 ± 0.6 


1.6 ± 0.2 


L S 


07 36 34.0 


+65 


35 


40.7 


0.7 ± 0.3 


0.3 ± 0.1 




07 36 34.1 


+65 


38 


54.7 


54.4 ± 3.9 


21.4 ± 1.6 


L 


07 36 35.6 


+65 


36 


08.4 


3.1 ± 0.5 


1.2 ± 0.2 




07 36 37.9 


+65 


38 


15.7 


2.0 ± 0.8 


0.8 ± 0.3 




07 36 37.9 


+65 


37 


37.2 


1.2 ± 0.4 


0.5 ± 0.2 




07 36 42.0 


+65 


36 


51.8 


52.6 ± 2.0 


20.7 ± 0.8 




07 36 43.0 


+65 


37 


09.7 


0.7 ± 0.3 


0.3 ± 0.1 




07 36 45.2 


+65 


35 


57.3 


3.1 ± 0.5 


1.2 ± 0.2 




07 36 45.8 


+65 


36 


40.5 


2.4 ± 0.5 


0.9 ± 0.2 




07 36 46.1 


+65 


36 


13.6 


-0.0 ± 0.1 


-0.0 ± 0.1 


L 


07 36 46.8 


+65 


35 


58.1 


1.0 ± 0.3 


0.4 ± 0.1 




07 36 47.5 


+65 


36 


19.0 


7.3 ± 0.8 


2.9 ± 0.3 


L S 


07 36 47.6 


+65 


36 


22.9 


5.5 ± 0.7 


2.1 ± 0.3 




07 36 50.1 


+65 


36 


03.7 


20.9 ± 1.3 


8.2 ± 0.5 


L 


07 36 52.0 


+65 


36 


41.0 


2.5 ± 0.5 


1.0 ± 0.2 




07 36 52.5 


+65 


36 


46.0 


0.9 ± 0.3 


0.4 ± 0.1 




07 36 53.8 


+65 


33 


32.1 


0.9 ± 0.2 


0.4 ± 0.1 




07 36 53.9 


+65 


35 


35.2 


0.5 ± 0.2 


0.2 ± 0.1 




07 36 55.4 


+65 


36 


08.3 


2.5 ± 0.4 


1.0 ± 0.1 




07 36 55.6 


+65 


35 


40.9 


436.8 ± 5.0 


172.3 ± 2.0 


L S 


07 36 56.1 


+65 


37 


15.6 


2.4 ± 0.5 


0.9 ± 0.2 




07 36 57.2 


+65 


36 


03.7 


2.8 ± 0.4 


1.1 ± 0.2 




07 37 00.7 


+65 


36 


06.2 


1.1 ± 0.3 


0.4 ± 0.1 




07 37 00.8 


+65 


34 


17.9 


14.2 ± 0.9 


5.6 ± 0.4 


L S 


07 37 01.8 


+65 


38 


27.2 


-0.1 ± 0.2 


-0.0 ± 0.1 




07 37 02.2 


+65 


34 


45.8 


0.4 ± 0.2 


0.2 ± 0.1 




07 37 02.4 


+65 


36 


01.5 


1.2 ± 0.3 


0.5 ± 0.1 




07 37 02.6 


+65 


37 


10.8 


2.4 ± 0.4 


0.9 ± 0.1 




07 37 02.9 


+65 


34 


37.9 


1.4 ± 0.3 


0.6 ± 0.1 




07 37 06.4 


+65 


34 


51.9 


0.7 ± 0.2 


0.3 ± 0.1 




07 37 06.9 


+65 


36 


21.7 


0.7 ± 0.2 


0.3 ± 0.1 




07 37 07.0 


+65 


36 


14.1 


1.1 ± 0.3 


0.5 ± 0.1 




07 37 07.1 


+65 


35 


56.6 


2.1 ± 0.4 


0.8 ± 0.1 




07 37 07.4 


-1-65 


34 


55.9 


21.1 ± 1.1 


8.3 ± 0.4 


L S 


07 37 09.2 


-1-65 


35 


44.2 


23.8 ± 1.2 


9.4 ± 0.5 


L 


07 37 09.4 


+65 


35 


01.4 


0.3 ± 0.2 


0.1 ± 0.1 




07 37 09.6 


+65 


33 


05.8 


3.0 ± 0.4 


1.2 ± 0.2 




07 37 09.9 


+65 


35 


46.9 


0.7 ± 0.2 


0.3 ± 0.1 




07 37 10.2 


+65 


33 


11.2 


3.5 ± 0.5 


1.4 ± 0.2 




07 37 11.3 


+65 


38 


14.2 


1.0 ± 0.3 


0.4 ± 0.1 




07 37 11.6 


+65 


33 


45.9 


51.7 ± 1.7 


20.4 ± 0.7 


S 


07 37 13.1 


+65 


35 


58.2 


3.0 ± 0.4 


1.2 ± 0.2 




07 37 14.9 


+65 


34 


29.0 


23.7 ± 1.2 


9.4 ± 0.5 


S 


07 37 16.1 


+65 


33 


28.9 


1.9 ± 0.3 


0.8 ± 0.1 




07 37 17.1 


+65 


35 


57.9 


42.4 ± 1.6 


16.7 ± 0.6 


L S 


07 37 17.6 


+65 


36 


23.4 


1.3 ± 0.3 


0.5 ± 0.1 




07 37 17.9 


+65 


37 


26.5 


1.2 ± 0.3 


0.5 ± 0.1 




07 37 18.0 


+65 


35 


09.4 


1.0 ± 0.3 


0.4 ± 0.1 





2.1. X-ray Source and Region Definitions 

The source-finding tool described by iTennantI ()2OO60 
was applied to all five datasets to search for discrete 
X-ray sources. The search was limited to the central 
6'x6' region and to events within the full Chandra en- 
ergy range 0.3—8.0 keV. Fifty-eight sources were detected 
in the merged dataset with a S/N above 2.8 and with a 
minimum of 5 cr above background uncertainty (corre- 
sponding to a detection limit of about 8—10 counts for 

^ http://cda.harvard.edu/chaser/ 



Note. — Units of right ascension are hours, minutes, and seconds, 
and units of declination are degrees, arcmintes, and arcseconds. 

Count rates for ASIC-S3, taking the 0.5 - 2.0 keV energy band. 
^ -Lx fof the 0.5 — 2.0 keV energy band. Lx is computed using the 
bremsstrahlung model obtained from Table |3] (see text). 

a typical on- axis source). This tally included all sources 
detected in the individual observations (provided they 
were within the FOV of the individual images) a nd al l 
those previously reported by iSchlegel fc Pannutil ()2003f ) 
within our 6'x6' field. The X-ray point sources and their 
properties are listed in Tabled 
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Since we are primarily interested in star-formation and 
the X-ray properties of recent and current star-forming 
events in this study, we also examined the X-ray emission 
from known SNRs and Hii regions. For this purpose, we 
first masked out the X-ray-detectcd point-like sources in 
the field of the merged image using a mask radius equal 
to 3 times the half-width of a Gaussian approximation 
to the instrumental point spread function (PSF). Larger 
radii were used, as necessary, for a few bright sources 
with strong emission in their PSF wings. Spectral infor- 
mation was then extracted from events falling within the 
area defi ned by optically identified SNRs from [SII]/Ha 
ratio by iMatonick et al.l ()1997[ ) using their listed loca- 
tions and source radii. One arcsecond ( ~ 16 pc ) radii 
are used for the SNRs whose sizes are unlisted. Local 
backgrounds were also extracted from surrounding an- 
nuli. There are 29 cataloged SNRs within our 6'x6' field. 
There are also three candidate radio SN Rs in our field 
(jPannuti et al.ll2007| ). One (denoted ^ bv iPannuti et al.l 
2007) is already included in the lis t of optically-id entified 
SNRs. Another (denoted TH2 bv IPannuti et al.l 2007) is 
a strong X-ray source, Lx = 2 x 10 '^^ erg s~^, as previ- 
ously noted bv IPannuti et al.l (|2007[ ). However, no radio 
spectral indices are available for any of these candidate 
radio SNRs. Therefore, we opt to omit TH2 from our 
list of candidate SNRs as it is anomolously bright if it is, 
in fact, a SNR. It remains classified as a discrete X-ray 
source in our tabulation. 

A similar proce dure was att e mpte d using the Hii re- 
gions defined by iSivan et aLl ()199Gf ) based on photo- 
graphic images but many of these cataloged objects 
did not correspond to Hii regions visible in a recent 
continuum-subtracted Ha CCD image. This is due 
mostly to crude modeling of the Hii structure which re- 
sulted in poor estimates of source positions and sizes. 
Therefore, we constructed our own catalog of Hii regions 
defined as circular approximations to the areas enclosed 
by surface brightness contours corresponding to a level of 
20 X the background in our continuum-subtracted Ha im- 
age. There are 58 Hii regions defined in this way within 
our 6' X 6' field includin g 5 of the 6 giant Hii regions iden- 
tified previously bv iDrissen et al.l (|1999i) . 

There are four point-like sources detected by our 
source-finding trial that spatially coincide with SNRs. 
Two more detected sources lie within a bright extended 
Hii region. Upon further examination, all these sources 
were determined to be steady thermal sources and were 
therefore assigned to the SNR (Hii) category rather than 
to the general category of point sources for purposes of 
our analysis. Conversely, the PSF wings of bright point 
sources overlap with the positions of two of the cataloged 
SNRs. These two SNRs were therefore excluded from our 
analysis because of this contamination. The three his- 
torical supernovae, SN 2004dj (a strong X-ray source), 
SN 2002kg and SN 1954J (both fainter than our detec- 
tion limit) and the X-ray-bright candidate radio SNR 
were also not considered further in this study. 

Finally, we also examined the properties of the X-ray 
emission remaining after the X-ray point sources and the 
optically-identified SNRs and Hii regions are excluded. 
This residual component is likely comprised of unresolved 
(faint) XRBs and other low-luminosity objects related to 
the stellar content of NGC 2403, extended diffuse hot gas 



TABLE 3 
Point sources fitting results 



fit Parameters Brcmsstrahlung 


Powerlaw 


A^H (102icm-2) 




■^■■^ -0.2 


kT (keV) 


4 7+0.6 

^•'-0.5 




r 




-I 07+0. 06 


Normalization (10^^) 


o.y_o 2 


^ 4+0-3 


Lx(0.5-2.0 keV)/1038 erg s"! 


7 7+0.3 
'•'-0.4 


7 7+0.5 
'■'-0.5 


Lx(2.0-10 keV)/1038 erg s"! 


16.7l?-» 






185/197 


211/197 



within the disk and halo of NGC 2403, and unresolved 
background (primarily AGN) and foreground (local dif- 
fuse Galactic) emission. Inspection of the residual emis- 
sion image shows a clear excess above the background in 
the central regions of NGC 2403. 

We examine the X-ray properties of these four source 
types in the following subsections. 

2.2. Discrete X-ray-Detected Sources 

iSchlegel fc PannuTil ()2003f ) have already presented the 
properties of the discrete source population detected in 
the first Chandra observation of NGC 2403. Here, our 
interest is only in the average spectral shape of the dis- 
crete sources and their luminosity distribution. The aver- 
age spectrum can be represented by co-adding the source 
spectra from only the first observation. This provides 
sufficient source counts for characterizing the spectrum 
yet avoids complications that arise when combining ob- 
servations made at different times during the Chandra 
mission. The spectra of the individual point sources (and 
corresponding backgrounds) were co-added and weighted 
ancillary response (ARF) and response matrix (RMF) 
files were generated using CIAO tool acisspec. The 
resulting spectrum was grouped to ensure at least 20 
counts in each spectral bin. The source spectrum was 
modelled using the XSPEC spectral-fitting utility (ver- 
sion 11.3. 2t). 

We applied both absorbed powerlaw and 
brcmsstrahlung models to the co-added discrete 
source spectrum. The fitting results are shown in Ta- 
ble [HI and the spectrum is shown in Figure [H Although 
we are mainly interested in the energy range 0.5—2.0 
keV, a more precise fit can be made by extending the 
spectral fits to at least 5 keV. Over this range, the 
best-fitting model is the thermal brcmsstrahlung model. 
This model is used to estimate the luminosities of the 
individual sources. For this purpose, separate ARF 
and RMF files for each source and for each observation 
were created and used in the fit. The hydrogen column 
density and brcmsstrahlung temperature were held 
fixed at the values from Table [3] during the fitting; 
only the model normalization was allowed to vary. The 
resulting 0.5—2.0 keV luminosities, averaged over the 
four observations, are listed in Tabled 

The flux estimates were used to check for variability of 
sources between and during individual Chandra observa- 
tions. This check helps to differentiate between diffuse 
emission associated with star-formation that has a rela- 
tively soft X-ray spectrum and no variability and emis- 
sion from XRBs that is spectrally hard and often variable 
over time. 
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Fig. 2. — Top: The co-added (stacked) spectrum of all 58 
discrete sources detected in the Chandra X-ray data. The co- 
added spectrum was constructed using only the first observation. 
Also shown are the best-fit model (an absorbed bremsstahlung) 
and the fit residuals. Bottom: The X-ray luminosity func- 
tion on the energy range of 0.5—2.0 keV (solid line) for the 58 
sources. The best-fit model was used to convert count rates to 
flux (3.1xl0~* ct s~^ corresponds to 10~^^ergs cm~^ s~^). The 
single powerlaw model fitted to the luminosity function for sources 
brighter than 2x10"^'' ergs cm~^ s^^ ~2xl0''^erg s~^ is shown 
as a solid curve and gives a slope of —0.54. The dashed curve rep- 
resents the expected contribution to the luminosity function from 
cosmic X-ray background sources as given by the Moretti ct al. 
(2003) analytic expression. 



If the flux from a source during one observation devi- 
ated from its average over the other available observa- 
tions by >3(T, then it was flagged as long-term variable 
and designated by the letter "L" in column 5 of Table [21 
(Because of differences in spacecraft roll and target aim- 
point, not all sources were within the field of view during 
all four observations.) Short-term variability during in- 
dividual observations was checked by binning the source 
light curve (event arrival time) for each observation us- 
ing 2 ks temporal bins. The shape of the light curve was 
compared to a constant using the statistic to check for 
variability. Short-term variability is also given in column 
5 of Table [2] by the letter "S". 

The cumulative X-ray luminosity function (XLF) for 
the 58 sources detected within the 6'x6' central re- 
gion of NGC 2403 is shown in Figure [2l Luminosi- 
ties for each source were computed from their aver- 
age flux over up to four observations, depending on 
availability, on the energy range 0.5—2.0 keV. Using 
only the 50 sources brighter than a flux limit flux of 



2.0x10 ergs cm ^ s ^ in the 0.5—2.0 keV energy 
range (corresponding to a luminosity at the distance 
of NGC 2403 of 2.4x10^^ erg s'^), a single power 
law gives a best-fitting slope of —0.54 ± 0.02, roughly 
consistent with the value of —0.59 ± 0.02 reported by 
iSchlegel fc Panniitil (|2003D . For reference, a flux of 
10~^^ ergs cm~^ in the 0.5—2.0 keV band corre- 
sponds to a flux of 2.1x10""'^^ ergs cm~^ s"^ in the 2.0- 
10.0 keV band assuming our best-fit co-added discrete 
source spectral shape. 

An estimate of the contribution from unrelated back- 
ground AGN to the detected source population in the 
field can be made by comparing the observed XLF to the 
analytical fi ts to the deep field cosmic X-ray background 
logTV-logS" fMorett i et "al]|200 3'). This fit is shown in Fig- 
ureO An estimated 21 ± 1 of the 50 sources, or 42±5%, 
are likely background AGNs. 

2.3. Supernova Remnants 

Most of the individual SNRs (and Hii regions discussed 
in the next subsection) are too faint to be detected as dis- 
crete X-ray sources. The possibility exists that many of 
the X-ray events detected within the spatial regions de- 
fined by the individual optical SNRs are unrelated "back- 
ground" events. Figure [3] shows the XLF for the SNRs. 
Also shown is a pseudo-random sampling of the back- 
ground defined using the same locations and sizes of the 
cataloged SNRs but with their RA coordinates reflected 
about the center of NGC 2403. (For both populations, 
only those regions with >1 background-subtracted net 
counts are shown.) This second XLF indicates that only 
the (X-ray) brightest 4 or 5 cataloged SNRs are truly 
sources of X-ray emission. Fitting a single power law 
model to the XLF of the eight SNRs with >4 source 
counts in the 0.5—2.0 keV energy range (corresponding 
to 1.1x10""'^^ ergs cm~^ s~^) gives a best-fitting slope of 
-0.89 ± 0.14. 

There are four optically-identified SNRs detected in 
X-rays via the source-finding tool. The X-ray-brightest 
SNR is also the brightest in Ha line emission. However, 
there is no overall correlation between X-ray and Ha lu- 
minosity in the SNR sample. Th is is consistent with the 
analysis of iPannuti et al.l (I2007D. Using the two S NRs 
whose diameters are listed in iMatonick eFall (fT99l . we 
also find the mean diameter of X-ray-detected SNRs is 
smaller (40 pc) than the mean diameter of X-ray non- 
detect ions (70 pc) again comfirming the IPannuti et al.l 
(|2007f) result. 

A co-added spectrum of the cataloged SNRs was cre- 
ated by merging the first three Chandra images (the 
fourth observation does not adequately cover the cen- 
tral 6'x6' region of interest). Twenty-four SNRs (of 27 
total) are imaged on S3 in the first three observations. 
The same averaging method used for the discrete sources 
(§ I2.2p was used here to build the co-added spectrum. 
Weighted ARF and RMF files generated from the sec- 
ond observation were used for the spectral fitting. We 
checked that this did not adversely affect our fitting re- 
sults by repeating the spectral analysis using ARFs and 
RMFs generated from the first and third observations. 

We used a surrounding annulus to define a local back- 
ground for the individual SNRs, after masking out the 
detected X-ray point sources. The total number of X- 
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Fig. 3. — Top: The co-added (stacked) spectrum of the 24 
optically-identified SNRs imaged on S3. The best-fitting model 
(an absorbed apec model with hydrogen column density, A^h, al- 
lowed to be a free parameter) is shown along with the fit resid- 
ual. The spectrum has been re-binned to have a width of 73 
eV for display purposes. Bottom: The luminosity function of 
the background-subtracted SNRs (solid line) and of the pseudo- 
random positions (dashed line, see text). The single powerlaw 
model fitted to the luminosity function of sources brighter than 
1.1x10"^® ergs cm~-^ ~2xl0^^erg is shown as a curve 
and gives a slope of —0.89. 



TABLE 4 

SNRs X-RAY FITTING RESULTS 



Fit Parameters 


apec 


apec 


apec 


A^H (lO^^cm-^) 
Te (MK) 
Normalization" 
Lx/lO^'^ erg s-l 

Lint ^Iq37 gj,g g-1 

C-stat/bins 


0.4 

o.i«- 

0.3 

107.4/108 


1.4 

3+°-l 
0.6 

105.9/108 


o g+:^. 3 
^■^-2.3 
n pr + 2.1 
o.0_1.4 

5+2-3 

Q O + 0.5 

0.8 

105.3/108 


Derived Parameters 


b 




V [/] (10"^ cm^) 


0.13 


0.13 


0.13 


ne (cm-3) 


0.04 


0.05 


0.07 


■P/fc [/"^^^] (10^ K cm-^) 


3.5 


4.1 


5.1 


(105 Mo) 
Eth (1053 erg) 


0.05 
0.09 


0.07 
0.1 


0.11 
0.1 


t. (Myr) 


48 


30 


16 



K : 



ray events in the co-added spectrum between 0.4 and 
2.0 keV is 153. This reduces to 94.8 ± 13.1 after back- 
ground subtraction. The XSPEC model apec, represent- 
ing thermal emission from an optically-thin coUisionally- 
ionized plasma, combined with a model phabs for an 
intervening absorption column were applied to the spec- 
trum. The parameters of the apec model are the plasma 
temperature, elemental abundances, and emission inte- 
gral. Elementa l abundances were f i xed to the solar val- 
ues as given in An ders fc Grevessd ([1989'). Three model 
fits were attempted as detailed in Table [H These dif- 
fered in the values assigned to the absorbing columns. 
In two of these models, the absorbing column was held 
fixed at either the Galactic value along the line of sight, 



4x10^° cm 



^, implying no local absorption within 
= 14 X 10^° cm~^ which is equivalent 



''volume filling factor scaling in [ ] in Column 1. 



NGC 2403, or iVn = 14 x lU'" cm 
to the mean value of the optical extinction derived for ac- 
tive star-forming regions (§[3]). The third model allows 
the absorbing column to be a free parameter in the fit. In 
this and in models of the Hii regions and of the residual 
X-ray emission to be discussed below, allowing the ab- 
sorbing column, plasma temperature, and normalization 
(emission integral) to vary in the fitting process often 
leads to large uncertainties in some of these parameter 
values. The reason for this is simply that the plasma 
temperature is quite low, of order a few 10^ K (a few 
0.1 keV), so that most of the emission is at low energies 
where the sensitivity to absorption is most acute. This 
introduces a degeneracy in which an equally acceptable 
fit is possible by increasing the column density, TVh, while 
simultaneously increasing the model normalization, usu- 
ally with a moderate decrease in the fitted temperature. 

The X-ray luminosities of the individual SNRs were 
computed by adopting the spectral model that best fits 
the stacked spectrum (Table lU column 4) . This model 
was applied to the distribution of events of each indi- 
vidual SNR using the best-fit values of and A^h, and 
allowing only the model normalization to vary in the fi- 
nal fitting. This is equivalent to scaling by the number 
of counts in the individual spectra. The X-ray count 
rates, corresponding luminosities, and Ha luminosities 
are listed in Table [5l 

A number of other physical parameters can be derived 
from the fitted values of the plasma temperature and 
emission integral provided some estimate of the volume 
occupied by the hot gas can be made. The total volume 
of the 24 SNRs is 1.3 x 10^^ cm'^ assuming a spheri- 
cal geometr y for each remn ant and adopting the radii 
reported in IMatonick et al.] (1997; the same radii used 
here to define the X-ray spectral source regions). The 
actual volume occupied by the hot gas is some fraction, 
/, of this total volume. Assuming the ion and electron 
number densities are equal, which is roughly true for a 
hydrogen-dominated hot plasma, then the emission in- 
tegral / neUndV ^ n^fV so that the electron density 
He oc {K/ fVy/^ where K is the spectral model nor- 
malization parameter which is itself proportional to the 
emission integral. 

From Ue and the flux- weighted mean temperature Te, 
we can estimate the hot gas pressure P/k = 2neTe, mass 
Mx = nefJ-rUpV where nip is the proton mass and /i = 1.4 
is the mass per proton, thermal energy Eth = SueVkTe 
and the cooling time tc = Eth/ L\jo\ of the X-ray emitting 
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TABLE 5 
SNRs 



RA Dec Count Rate " Lx " Lua 

(J2000) (J2000) 10-5 cts lO^Serg s'^ lO^^erg s-^ 

07 36 42.8 +65 34 52.7 -0.8 ± 1.0 -0.3 ± 0.4 0.3 

07 36 45.7 +65 36 35.7 1.1 ± 1.5 0.4 ± 0.6 1.7 

07 36 45.7 +65 36 40.7 25.8 ± 4.8 10.0 ± 1.8 17.1 

07 36 49.0 +65 34 31.0 3.3 ± 3.8 1.3 ± 1.5 2.4 

07 36 52.0 +65 33 41.0 -0.5 ± 0.9 -0.2 ± 0.3 0.5 

07 36 52.7 +65 35 50.5 -0.8 ± 1.0 -0.3 ± 0.4 4.8 

07 36 53.3 +65 36 00.3 3.8 ± 2.0 1.5 ± 0.8 1.0 

07 36 53.5 +65 33 42.0 2.4 ± 2.1 0.9 ± 0.8 5.5 

07 36 53.6 +65 35 11.6 -0.5 ± 0.9 -0.2 ± 0.3 2.9 

07 36 56.2 +65 34 06.0 5.1 ± 2.7 2.0 ± 1.1 1.8 

07 36 57.0 +65 36 04.8 -0.5 ± 0.9 -0.2 ± 0.3 0.3 

07 37 01.6 +65 34 13.0 7.8 ± 3.4 3.0 ± 1.3 2.9 

07 37 02.0 +65 34 37.0 0.0 ± 0.9 0.0 ± 0.4 5.5 

07 37 02.0 +65 33 42.0 0.3 ± 0.9 0.1 ± 0.3 0.4 

07 37 02.2 +65 37 20.7 1.0 ± 1.2 0.4 ± 0.5 2.3 

07 37 02.2 +65 36 02.0 3.7 ± 1.8 1.5 ± 0.7 2.9 

07 37 03.0 +65 33 45.0 0.0 ± 0.9 0.0 ± 0.4 0.3 

07 37 03.0 +65 34 38.0 10.1 ± 2.9 3.9 ± 1.1 1.0 

07 37 06.0 +65 36 03.2 -2.4 ± 1.8 -0.9 ± 0.7 2.0 

07 37 06.1 +65 36 10.3 -0.8 ± 1.7 -0.3 ± 0.7 2.2 

07 37 07.1 +65 37 10.6 -1.3 ± 1.0 -0.5 ± 0.4 1.2 

07 37 10.6 +65 33 10.0 2.8 ± 1.9 1.1 ± 0.7 0.3 

07 37 12.5 +65 33 46.0 0.5 ± 0.8 0.2 ± 0.3 2.0 

07 37 16.1 +65 33 29.0 13.2 ± 3.3 5.1 ± 1.3 0.8 

Units of right ascension are hours, minutes, and seconds, and units 

of declination are degrees, arcmintes, and arcseconds. 
" Count rates for ACIS-S3, talking the 0.5 — 2.0 keV energy band. 

Luminosity in the 0.5 — 2.0 keV energy band. Lx is computed 
using the one temperature model with variable model as listed 
in Table |4] (see text) 



plasma as listed in Table [H Here, Lboi is estimated by 
integrating the spectral model from 0.02 keV to 20 keV 
within XSPEC. For plasma temperature typical of SNRs 
(and of other thermal sources analyzed in this paper), 
Lboi is 5 times Lx, where Lx is the observed X-ray 
luminosity quoted on the standard 0.5—2.0 keV range. 
These derived quantities are all functions of the volume 
filling factor, / as indicated in the table. 

The total X-ray luminosity from the 24 SNRs is only 
3x10'^^ erg s~^. Inspection of Figure [3] shows that most 
of this emission comes from the few brightest sources 
with the most luminous being ~lxl0^^ erg s~^. The 
thermal energy content of the SNRs, ~10^^/^/^ erg, cor- 
responds to up to ~40% (in the limit / = 1) of the initial 
kinetic energy released by the supernova explosions, as- 
suming 10^^ erg per event. The total mass of hot gas is 
substantial. It corresponds to an average mass per SNR 
of Mx ^ 300 ± 100/1/2 M© (uncertainties given by the 
range of values deduced from the three X-ray spectral 
models). This can be compared to the amount of ma- 
terial swept up by the SNRs, ^ 5000 nu per SNR, 
assuming a uniform density uh cm^'^, over the SNR vol- 
ume. 

2.4. Hll Regions 

We performed an analysis of the 48 Hii regions, im- 
aged on S3 in the first three observations, analogous to 
that of the SNRs. Figure H shows the XLF for the Hii 
regions and the pseudo-random sampling defined using 
RA coordinates reflected about the center of NGC 2403. 
Again, only the most X-ray-luminous Hii regions are 
significant detections. Similar to the SNRs, we also 
fitted a single power law model to the XLF of the 7 



brightest Hii regions, (>4 counts and corresponding to 
1.1x10^1^ ergs cm~2 s"^). The best- fitting slope is 
—0.47 ± 0.13, somewhat shallower than for the SNRs. 

The co-added X-ray spectrum is shown in Figure |4l 
Similar to the SNR spectrum, we also applied thermal 
models to the Hii region spectrum. There are sufficient 
source counts (858 cts, 422±32 after the background 
subtraction) in the Hii region spectrum to allow 2- 
temperature model fits in addition to the 1-temperature 
models that were applied to the SNR spectrum. The fit- 
ting results and derived values are listed in Table [S) We 
assumed pressure equilibrium between the two temper- 
ature components and identical filling factors to derive 
the values of ne, i?t/i, and Mx- 

Note that the best-fitting 1-temperature model with 
variable column density results in a very high model nor- 
malization which leads to unreasonable values of the de- 
rived parameters. However, this does not seem to be the 
issue with the 2-temperature model. This model (col- 
umn 6 in Table (5]) results in A^h ~ A^h (though with 
large uncertainty) and hence with values of the remain- 
ing parameters similar to those of the fixed A^h model. 

The X-ray luminosities of the individual Hii regions 
were computed in the same manner as for the SNRs; us- 
ing the best-fitting two temperature spectral model with 
variable A'h in this case. The X-ray count rates, corre- 
sponding luminosities, and Ha luminosities are listed in 
Table [3 

There are twice as many Hii regions in our sample 
as there are SNRs in our SNR sample. The X-ray 
temperatures of the Hii regions are about 40% lower 
compared to the SNRs. The total X-ray luminosity is 
^1.5x10'^^ erg s"^ for the Hii regions which is only a 



8 



— — \\ 












1 



0.5 1 2 

channel energy (keV) 

T 1 , 1 I 1 I ' ' ' 




0.1 1 10 100 

(10"^'' erg"' cm"^ s"') 



Fig. 4. — Top: The stacked spectrum of the identified Hll regions. 
The best-fitting model (an absorbed two-temperature apec model 
with hydrogen column density, Nn, as a free parameter) is shown 
along with the fit residuals. The spectrum has been re-binned to 
have a width of 73 eV for display purposes. Bottom: The luminos- 
ity function of the background-subtracted Hll regions (solid line) 
and of the pseudo-random positions (dashed line, see text). The 
single powerlaw model fitted to the luminosity function of sources 
brighter than 1.1x10"^^ ergs cm~^ s~^ ~2xl0^^erg s~^ is shown 
as a curve and gives a slope of —0.47. 



factor of 3 higher, on average, for individual Hii regions 
compared to the SNRs. Similarly, the average hot gas 
mass is r-^1000 Mq, per Hii region, based on the best- 
fitting 2-temperature spectral model. This is about a 
factor of 3 larger than the hot gas mass in individual 
SNRs. 

2.5. Residual Emission 

The left panel of Figure displays a smoothed image of 
the soft X-ray emission from the underlying unresolved 
component. This residual emission is defined as the net 
emission left after masking out all detected point sources 
(including SN 2004dj and TH2), Hii regions and SNRs. 
For spectral fitting, an X-ray map of this component was 
created using only the first Chandra observation. This is 
a compromise in that a deeper image could be made using 
an image merged from two or more individual observa- 
tions but the area of the region of overlap sampled by 
these multiple observations is too small to define both 
a sizable source region and a surrounding background 
region that are both wholly contained within the over- 
lap region (see Figure [l}. The first Chandra observation 
placed the center of NGC 2403 near the S3 aimpoint and 



thus is the most useful for our purposes. The source 
region is defined as a disk of radius 2.'5 (2.3 kpc) cen- 
tered on NGC 2403. The background is composed of 
three non-contiguous regions of low surface brightness 
as shown in Figure [1] To construct the smoothed im- 
age depicted in Figure [5l all excluded regions were filled 
with a local background level using the Poisson method 
of the CIAO tool dmf ilth^, then divided by an exposure 
map evaluated at an energy of 0.5 keV. This exposure- 
corrected image was then smoothed using the CIAO tool 
aconvolve using a Gaussian function with a width of 10 
pixels (~ 5"). 

There are several nearly point-like regions that remain 
in the residual emission image. These may be faint point 
sources that were hidden in the PSF wings of brighter 
nearby sources or extended non-spherical features in Hii 
regions not covered by our (circular) masks. We exam- 
ined the five brightest of these regions and find they con- 
tribute only 4% of the counts in the residual emission 
spectrum. Therefore, these regions were left in as part 
of the residual emission. 

In general, the residual emission fills a broad region 
of the central part of NGC 2403 extending from north- 
west of the nucleus (located at the center of the image) 
to southeast of the nucleus. This region is not a par- 
ticularly strong source of emission at UV, mid-IR, nor 
Ha wavelengths. For example, a grayscale UV image 
is shown in the right panel of Figure \5\ with the X-ray 
contours overlaid. This image suggests that much of the 
residual X-ray emission is centrally located relative to 
the UV-bright zones that contain stars ^ 100 Myrs old. 

Several models were fitted to the spectrum of this resid- 
ual X-ray emission. Area-weighted ARFs and RMFs 
were created for this purpose using the CIAO script 
specextract. The spectral fitting results are listed in 
Table [i 

The three thermal models with the absorption column 
density frozen during the fitting all give low values for 
the model normalization, electron density, pressure, hot 
gas mass, and thermal energy. The best-fitting of these 
models is the 2-tempcrature model with temperatures 
of 2.3 and 8.6 MK with most of the flux in the cooler 
component (Tabled column 5), shown in Figure [G] Fix- 
ing the absorption column density to the Galactic value 
implies that the residual emission lies above the disk of 
NGC 2403 where there is little or no overlying cold gas. 
This is a reasonable assumption if the X-ray-emitting gas 
was ejected from the galaxy through galactic-scale winds 
or smaller fountains or chimneys. Fixing the absorption 
column density to the equivalent mean value of the op- 
tical extinction toward the active star-forming regions 
in NGC 2403, Nu = 1.4 x 10^^ cm'^, corresponds to a 
modest layer of overlying cold gas. 

The models with the absorption column density free 
to vary during the fit result in statistically improved 
fits with absorbing columns higher than found for any 
of the other types of X-ray sources considered in the 
previous subsections. Such high columns imply that 
the hot gas is located behind a layer of neutral gas 
such as would be the case if the hot gas were con- 
fined to the disk. However, the observed Hi col- 
umn densities through the central disk of NGC 2403 

^ See jhttp: / / cxc.harvard.edu/ciao/tlireads/diffuse_emission/| 
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TABLE 6 

Hll REGIONS X-RAY FITTING RESULTS 



Fit Parameters 


apec 


apec 


apec 


apec+apec 


apec+apec 


iVn (102lcm-2) 
Te, (MK) 

Normalization Ki " 
Te, (MK) 


0.4 

^ =+0.4 

^•'-0.1 


1.4 

Q+0.6 

1 ^+0.3 


1 1+'^''' 
^•^-0.2 

91 t;q+33590 


0.4 


Q 4+2.3 
S.8t\i 


Normalization K2 
Lx/lO^'^ erg s"! 


1 9+0.2 
^•^-0.2 


1 9+0.2 
^■^-0.3 


1 n + 0.5 

-^•■^-1.0 


1 4+0-2 


3+°'2 

1 4+°l 


L™*/10^^ erg s-i 


1.5 


2.7 


365 


1.7 


1.7 


C-stat/bin 


154.0/108 


151.9/108 


129.4/108 


120.3/108 


120.3/108 


Derived Parameters'" 


y [/] (10"^ cm^) 


2.2 


2.2 


2.2 


2.2 


2.2 


ne (cm-3) 


0.02 


0.03 


1.1 


0.07,0.01 


0.07,0.01 


P/k (10^5 K cm"^) 


1.4 


1.9 


40 


2.5 


2.5 


(iqS M0) 


0.5 


0.8 


28 


0.5 


0.5 


Eth (10^3 erg) 


0.6 


0.9 


18 


1.1 


1.1 


[/+1/2] (Myr) 


55 


32 


0.1 


32 


32 



K = (10~y/47rD^) JneTipdV 
'volume filling factor scaling in [ ] in Column 1. 




Fig. 5.— The residual X-ray emission from NGC 2403. Left: The soft X-ray image (0.4-2.0 keV) of the central 6'x6' region. All 
four observations have been co-added and exposure-corrected and all discrete sources, optically-identified SNRs and Hll regions have been 
removed. The image has been aggressively smoothed with the CIAO tool aconvolve. The contours indicate 10, 20, 30, and 40(t above the 
background level. The gray scale indicates intensity levels in units of photon cm~^ s~^ pixel"'^. One pixel corresponds to 0"492. Right: 
The X-ray residual emission contours overlaid on the GALEX FUV image. The gray scale of the image indicates intensity levels in units 
of ct s~^ pixel~2. One pixel corresponds to 1"5. 



range from only '^l O^^ cni~^ ( Thornley fc Wilson |[r995fl 
to ^3 x10^^ cni-2 (jSchaap et^ 2000; Fratcrnai ret al.l 
l2002bf l which is considerably less than implied by these 
X-ray model fits, suggesting the models are unphysi- 
cal. Furthermore, inspection of the images of this hot 
gas component and comparison to images at other wave- 
lengths suggest that the hot gas is not emitted from re- 
gions with high levels of overlying cold gas. 

Therefore, the models with the absorption column 
fixed are the most realistic models. They imply (4 — 
7) X 10^ M0 of hot gas is present in the central regions of 
NGC 2403, assuming a 2.3 kpc radius spherical emission 
volume. This hot gas is rather tenuous, Ue ~ 0.003 to 
0.004 cni^-^, with, formally, a long cooling time of 200 
— 340 Myr. We note that the actual X-ray emitting re- 



gion is slightly smaller than the 2.3 kpc radius spectrum 
extracting region, so that the assumed volume may be 
larger than the volume containing the hot gas. If we as- 
sume instead that the hot gas is confined to a disk with 
thickness 200 pc and 2.3 kpc radius, then the inferred 
density increases by about a factor of 4 in compensation 
for the smaller volume for a given emission integral. This 
reduces the cooling time to 50 — 90 Myrs. The true emis- 
sion volume is probably somewhere between our sphere 
and disk estimates. Table [8] contains derived properties 
of the residual gas for both these geometries. 

An estimate of the contribution to this residual X-ray 
emission from unresolved sources fainter than, but oth- 
erwise similar to, the discrete source population ft j2.2p 
can be made by extrapolating the fitted point-source lu- 
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TABLE 7 

Hll REGIONS 



R.A. 


Decl 


Count Rate " * Lhci 


J2000 


J2000 


10-^ cts s lO^^'erg s-^ lO^^erg s"! 



07 36 35.7 


+65 


35 


09.0 


3.8 ± 


2.9 


1.6 


± 


1.3 


16.4 


07 36 38.0 


+65 


36 


09.5 


1.6 ± 


1.5 


0.7 


± 


0.6 


2.3 


07 36 38.4 


+65 


36 


23.7 


3.0 ± 


2.3 


1.3 


± 


1.0 


5.7 


07 36 39.6 


+65 


37 


11.5 


0.3 ± 


0.9 


0.1 


± 


0.4 


2.2 


07 36 41.7 


+65 


38 


06.7 


-0.2 ± 


1.4 


-0.1 


± 


0.6 


31.6 


07 36 44.6 


+65 


35 


07.7 


3.8 ± 


2.4 


1.6 


± 


1.0 


6.1 


07 36 45.6 


+65 


36 


59.9 


9.8 ± 


7.3 


4.2 


± 


3.1 


147.7 


07 36 45.7 


+65 


36 


37.4 


0.8 ± 


1.6 


0.4 


± 


0.7 


0.9 


07 36 45.9 


+65 


36 


40.6 


-0.3 ± 


0.8 


-0.1 


± 


0.3 


0.1 


07 36 45.9 


+65 


37 


14.6 


-0.8 ± 


1.9 


-0.3 


± 


0.8 


3.5 


07 36 46.7 


+65 


36 


36.9 


1.4 ± 


2.9 


0.6 


± 


1.2 


5.4 


07 36 46.9 


+65 


37 


05.3 


-1.7 ± 


2.5 


-0.7 


± 


1.1 


3.6 


07 36 47.1 


+65 


35 


42.5 


-1.9 ± 


1.8 


-0.8 


± 


0.8 


6.6 


07 36 47.2 


+65 


36 


44.8 


-0.8 ± 


1.7 


-0.3 


± 


0.7 


2.3 


07 36 47.4 


+65 


37 


00.1 


-0.4 ± 


1.9 


-0.2 




0.8 


1.2 


07 36 47.6 


+65 


35 


53.5 


1.1 ± 


1.8 


0.5 


± 


0.8 


2.5 


07 36 47.9 


+65 


34 


34.8 


0.0 ± 


0.9 


0.0 


± 


0.4 


2.2 


07 36 48.1 


+65 


33 


26.3 


0.0 ± 


1.6 


0.0 


± 


0.7 


7.0 


07 36 49.2 


+65 


36 


33.5 


-1.6 ± 


2.2 


-0.7 


± 


0.9 


8.5 


07 36 49.3 


+65 


36 


52.9 


-3.8 ± 


1.4 


-1.6 




0.6 


14.3 


07 36 50.0 


+65 


35 


11.2 


-0.5 ± 


1.4 


-0.2 




0.6 


4.8 


07 36 50.4 


+65 


35 


28.6 


-1.6 ± 


1.7 


-0.7 


± 


0.7 


2.8 


07 36 50.9 


+65 


34 


57.7 


-2.4 ± 


1.2 


-1.0 




0.5 


6.1 


07 36 52.3 


+65 


36 


47.3 


74.1 ± 10.3 


31.9 




4.4 


140.7 


07 36 52.9 


+65 


36 


20.7 


-2.1 ± 


1.2 


-0.9 


± 


0.5 


3.5 


07 36 53.2 


+65 


36 


34.7 


-0.8 ± 


1.0 


-0.3 




0.4 


0.8 


07 36 53.2 


+65 


34 


14.4 


-1.1 ± 


1.0 


-0.5 


± 


0.4 


3.8 


07 36 53.4 


+65 


35 


28.4 


-0.3 ± 


1.3 


-0.1 


± 


0.6 


1.2 


07 36 54.5 


+65 


35 


47.8 


-1.6 ± 


1.7 


-0.7 




0.7 


4.1 


07 36 55.6 


+65 


34 


11.0 


-0.3 ± 


1.9 


-0.1 


± 


0.8 


1.9 


07 36 57.8 


+65 


37 


24.5 


2.7 ± 


4.4 


1.1 




1.9 


67.1 


07 36 59.4 


+65 


35 


09.1 


0.8 ± 


1.9 


0.3 




0.8 


1.3 


07 36 60.0 


+65 


35 


06.8 


-2.4 ± 


1.4 


-1.1 


± 


0.6 


2.1 


07 37 00.3 


+65 


35 


11.0 


1.1 ± 


2.6 


0.5 


± 


1.1 


5.8 


07 37 00.8 


+65 


35 


25.2 


11.5 ± 


4.3 


5.0 


± 


1.9 


15.7 


07 37 03.7 


+65 


36 


19.9 


-0.8 ± 


0.9 


-0.3 


± 


0.4 


2.9 


07 37 05.8 


+65 


35 


52.3 


1.1 ± 


1.5 


0.5 


± 


0.7 


3.9 


07 37 07.0 


+65 


36 


41.1 


206.1 ± 16.6 


88.8 


± 


7.2 


448.2 


07 37 10.6 


+65 


35 


42.2 


-0.2 ± 


1.9 


-0.1 


± 


0.8 


6.3 


07 37 11.1 


+65 


34 


36.7 


-0.8 ± 


1.0 


-0.3 




0.4 


3.0 


07 37 11.6 


+65 


33 


56.2 


1.8 ± 


1.5 


0.8 




0.6 


1.5 


07 37 11.7 


+65 


33 


59.5 


0.9 ± 


1.6 


0.4 


± 


0.7 


2.4 


07 37 11.8 


+65 


33 


49.0 


-1.6 ± 


2.0 


-0.7 




0.8 


5.3 


07 37 14.0 


+65 


36 


10.0 


-1.9 ± 


1.1 


-0.8 




0.5 


1.7 


07 37 16.4 


+65 


34 


03.5 


1.3 ± 


2.4 


0.6 


± 


1.0 


5.3 


07 37 18.6 


+65 


35 


41.7 


1.2 ± 


2.5 


0.5 


± 


1.1 


8.2 


07 37 18.7 


+65 


33 


50.6 


8.3 ± 


3.7 


3.6 


± 


1.6 


36.2 



Note. — Units of right ascension are hours, minutes, and seconds, 
and units of decUnation are degrees, arcmintes, and arcseconds. 

Count rates for ASIC-S3, taking the 0.5 - 2.0 keV energy band. 
^ Luminosity in the 0.5 — 2.0 keV energy band. Lx is computed using 
the two temperature model with variable Nn model listed in Table |6] 



minosity function slope (Figure [2]) to lower luminosities. 
For this, we scale the XLF to include only those sources 
within the 2'.5 residual emission region. Using the same 
absorbed bremsstrahlung model as before 12. 2p results 
in an estimated unresolved point source contribution of 
2x10'^^ erg in the 0.5 to 2.0 keV energy range or 
about a 10-15% contribution to the diffuse emission. 
As an alternative estimate, we added a bremsstrahlung 
model component to the hot gas spectral model to esti- 
mate the contribution from unresolved sources. Adding 
a bremsstrahlung component does not improve statisti- 
cal significance of the fit but can contribute up to ^20% 
of the flux in the 0.5—2.0 keV band. The contributions 
from SNRs and Hii regions to the residual X-ray emis- 



sion can also be estimated from their XLFs (Figures [3l 
and[3]). These estimates are 9x10^^ erg s~^ for the SNRs 
and 1 X 10^^ for the Hii regions. The estimated contribu- 
tion from SNRs is relatively high because the XLF slope 
is steep for this source population. In any event, faint 
SNRs and Hii regions can contribute only a small frac- 
tion to the total X-ray luminosity of the residual emission 
which is l.GxlO'^^ erg s^^ or higher. 

2.6. Discussion of Preliminary Analysis of X-ray 
Observations 

The bulk X-ray properties of the SNRs, Hii regions, 
and diffuse (residual) emission are summarized in Ta- 
bles m [6l and [H respectively. They show that the hot 
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TABLE 8 

Residual emission X-ray fitting results 



Fit Parameters 


apec 


apec 


apec 


apec+apec 


apec+apec 


Nn (102icm-2) 

(MK) 
Normalization Ki " 
Te, (MK) 


0.4 
o 1+0.3 
'-'■^-0.2 

^•^-1.0 


1.4 

2.8ig:? 

17.2t\i 


5.6±i;'^ 

1 7+0.2 
-^•'-0.3 


0.4 

n q+0.4 
^■•^-0.3 
7 7+1.3 
'■'-1.2 

8.6tM 


9+3. a 

1 7+0.1 
^■'-1.7 
g49+2639 
"^^-432 
45 5+697 


Normalization K2 
Lx/lO^* erg s"! 
L™*/103* erg s-i 


1 O+0.2 
-^■'-'-0.2 

1.6 


-I Q+0.2 
^•'-'-0.2 

3.1 


1.5l?:? 
70 


9 c+o.e 

■^■^-0.7 
1 5+0.3 
^^■^-0.3 

1.9 


Q 9+13.8 
^•^-5.4 
1 7+0.4 
^- ' -0.9 

56.9 


C-stat/bin 


140.7/108 


133.9/108 


106.6/108 


105.9/108 


97.1/108 


Derived Parameters'" 




1535 


1535 


1535 


1535 


1535 


Tie (cm-3) 


0.003 


0.004 


0.026 


0.006,0.002 


0.08,0.003 


P/k [/"^/^] (10^ K cm-3) 


0.2 


0.2 


0.8 


0.3 


2.5 


[/+1/2] (105 Mo) 


47 


66 


460 


43 


167 


Eth (10^3 erg) 
(Myr) 


51 
311 


67 
203 


270 
7.3 


90 
342 


830 
30 


[/I (10'' cin') 


97 


97 


97 


97 


97 


ne (cm-3) 


0.01 


0.01 


0.10 


0.023,0.006 


0.31,0.01 


P/fc (105 K cm-3) 


0.7 


0.8 


3.3 


1.1 


10 


(105 Mq) 
(Myr) 


11 
12 

78 


17 
17 

51 


114 
66 
1.8 


11 
23 
86 


42 
211 
7.5 



For two tempareture models, the average is estimated as Mx / (fJ^rripV) and used in text. 

K = (10-9 /47rD2) f UeUpdV 
''Volume filling factor scaling in [ ] in Column 1. 
■^A 2.3 kpc radius spherical emission volume. 
"^A disk with the thickness 200pc and 2.3 kpc radius. 



SNR compared to about 1000 Mq per Hii region. Of 
course, Hii regions are composed of many stars a. n d SNs 
whereas the SNRs identified by iMatonick at all ()1997l ) 
are morphologically consistent with single events. Evi- 
dently, SNe occurring within the relatively tenuous and 
hot interior of Hii regions are much less efficient at pro- 
ducing hot gas than are isolated SNe. The ratio of the 
hot gas mass to the swept up mass for SNRs is there- 
fore ~0. 06— 0.10 and for the Hii regions it is only '^0.03 
for an assumed ambient density of 1 H cm^^. This indi- 
cates, again, that shock heating and thermal evaporation 
from the swept-up shells of young isolated SNRs is more 
efficient than in Hii regions. 

The diffuse (residual) emission has a density and pres- 
sure lower by about a factor of 10 (assuming the spherical 
volume) compared to SNRs and Hii regions but a tem- 
perature that is comparable to both these sources. If 
we assume a disk geometry, then the density is closer 
to the density in SNRs and Hii regions. There is an- 
other parameter to consider. The volume filling factor, 
/, may be lower for this residual emission compared to 
the SNRs and Hii regions because, by definition, this 
residual shares the galaxy disk with the much colder 
neutral and molecular gas and dust. Assuming / ~ 1 
for the Hii regions and pressure equilibrium suggests a 
filling factor / ~ (0.2 — 0.3) for the residual emission 
with correspondingly higher density and shorter cooling 
time. Alternatively, we may assume pressure equilib- 
rium between the hot residual gas component and the 
colder gas and dust. This cold gas resides primarily in 
clouds of typical temperatures of T ~ 100 K and den- 
sities of 60—100 cm^"^. This corresponds to a typical 
pressure P/k ~ lO'* cm^'^ K which would imply, at the 
observed temperature of (2 — 3) x 10^ K, a hot gas density 



> 

0) 




channel energy (keV) 



Fig. 6. — The X-ray spectrum of the residual emission from the 
central 6' X 6' region of NGC 2403. Data from only the first observa- 
tion are included. The best-fitting model, a two-temperature apec 
model with hydrogen column density fixed to the galactic value, 
and the fit residuals are shown. The spectrum has been re-binned 
to have a width of 73 eV for display purposes. 



gas densities in SNRs and Hii regions are comparable 
but that the SNR temperature is somewhat higher re- 
sulting in higher pressures, P oc nTx, by about a factor 
of 2. This is not suprising since the SNR, which mus t 
have ages of order 10^ yr or less (jMatonick et al.|[l997D . 
are expected to be further from pressure equilibrium with 
their surroundings than the fully-developed Hii regions 
which have ages of 10^ to >10^ yr. 

Another difference between the SNRs and Hii regions 
is the amount of hot gas estimated from the X-ray data. 
There is, on average, about 300±100 Mq of hot gas per 
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of ^0.003—0.005 cm ^ cons istent with the derived value 
(Table[8]) if / - 1. We note ()Fraternah et al.ll2002H) that 
the mean Hi density in the central regions of NGC 2403 
is about 0.4 cm~'^. Thus its filling factor must be only 
around 5% for a cloud density of 80 cm~'^ . A similar 
argum ent can be made for the H2 content (jSheth et al.l 
12001 ). Thus the filling factors for the cold gas may be 
quite small so that the bulk of the volume in the cen- 
tral regions of NGC 2403 may be filled with hot gas at 
high filling factor (/ ^ 1). This, in turn, would imply a 
considerable overpressure in the Hii regions which is not 
wholly unreasonable. 

In any case, the mass of hot gas in this residual com- 
ponent, Mx ^ 7 X 10^ Mq, is a small fraction of the 
mass contained in Hi (5x10^ M0) and comparable to 
the molecular (8x10^ Mq) gas mass. 

We cannot tell with certainty where this X-ray emit- 
ting residual gas is located. There is not a strong spatial 
correlation between this gas and star-formation indica- 
tors including Ha and UV emission. Combining this fact 
with an estimated cooling time of up to 340 Myr (de- 
pending on the geometry assumed and the filling factor) 
allows for several possibilities. The residual gas may have 
been created in situ during star-formation activity that 
ended up to 100 Myr ago (in regions currently lacking 
strong UV emission) or it may have been produced more 
recently but has escaped from its place of origin either 
into the halo or into low-density regions of the disk. 

At this point, we only have crude estimates of the ages 
of the various star-forming regions that are present in the 
central regions of NGC 2403. The ages of Hii regions are 
<10 Myr, equivalent to the lifetimes of the least-massive 
LyC-producing stars. The ages of UV-emitting regions 
can be as old as 100 Myr. What is needed is a better 
estimate of the age of individual star- forming regions. 
From this, we can better estimate how X-ray properties 
evolve through time and so better constrain the origin 
of the residual X-ray emission - whether it traces recent 
star-formation activity or is a relic of past activity. This 
is the subject of the next section. 

3. INDIVIDUAL STAR-FORMING REGIONS AND THEIR 
ANALYSIS 

We now turn to the analysis of individual young star- 
forming regions that have well-defined ages and masses. 
We choose young regions because they are more likely 
sources of X-ray emission from hot gas than are older 
clusters; and have more likely retained their identity as 
a coherent (though not necessarily bound) and coeval 
system under the prevailing tidal forces. Theoretically, 
X-ray emission from stellar winds and hot gas associ- 
ated with supernovae (SNe) occurs only within the first 
^40 Myr in the lifetime of a star cluster corresponding 
to the longest-lived stars that become core collapse SNe, 
roughly the 8—10 Mq stars. Star clusters in this age 
range emit strongly in UV light and (during the first 
'^10 Myr) in Ha line emission. However, these star- 
formation indicators are sensitive to extinction by overly- 
ing dust. Dust heated by UV and optical radiation emits 
this energy in the mid-IR and longer wavelength bands. 
Therefore, we define our regions as those bright in both 
UV radiation from massive young stars and in mid-IR 
radiation from reprocessing by dust. In this way, we can 
better estimate the extinction along the line of sight us- 



ing the mid-IR flux and thus correct the UV emission to 
obtain a better estimate of the age and mass of the star 
clusters. 

Star-formation tends to propagate spatially by com- 
pressing surrounding colder regions of the interstellar 
medium through the action of massive star winds and 
supernovae. Thus, small isolated regions are more likely 
to have a single characteristic age than are large ex- 
tended regions. Therefore, we apply our source-finding 
algorithm to identify isolated star-forming regions quan- 
titatively rather than to rely on a pre-determined ob- 
ject size (such as a fixed aperture or spatial grid, for in- 
sta nce, which is kn own to sample a range of stellar ages; 
cf.. iCalzetti et al.l 2005). The drawback to our selection 
method is that the flux from individual regions may be 
weak at wavelengths other than the UV and mid-IR we 
use to select them. 



3.1. Complementary Observations 

We use GALEX UV and Spitzer mid-IR measurements 
to identify young star- forming regions. We also use Ha, a 
traditional star-formation tracer, in our analysis to con- 
firm a young age for the regions. 

GALEX observed NGC 2403 on December 5, 2003 
Tilenum 5087) as par t of the Nearby Galaxy Survey 
Gil de Paz et al.l l2007f) . Corrected intensity maps at 
FUV (1529 A central wavelength) and NUV (2321 A) 
were obtained from the GALEX archive''. The images 
have a pixel scale of 1.5" and a spatial resolution of 
about 5". We applied the standard flux calibrations of 
1.40x10^15 ergs cm'^ s'^ cts'^ A'^ to the FUV and 
2.06 X 10-^^ ergs cm"^ g-icts-i A^^ to the NUV im- 
ages. 

NGC 2403 was observed with Spitzer IRAC and MIPS 
as a part of the Spitzer In frared Nearby Galaxies Survey 
(SINGS) legacy program (|Kennicutt et al.l[2003h . MIPS 
images at 24 /im were taken on 2004 October 13 and 
16 (key 5549568 and 5549824, respectively). The IRAC 
3.6 fiui, and 4.5 /im images were taken on 2004 October 8 
(key 5505792) and 12 (key 5505536). Final mosaiced im- 
ages provided by the SINGS program were used here. De- 
tails of their data reduction and calibration are reported 
in the SINGS data delivery paper®. The 24 /im image 
has an '^5" resolution which is comparable to GALEX. 
The final mosaiced image has a pixel scale of 1.5". Both 
the IRAC images have a spatial resolution of ^2" and a 
pixel scale of 0.75". 

As part of the SINGS project, NGC 2403 was ob- 
served with the Kitt Peak National Observatory 2.1m 
telescope in Ha bands on 2001 November 8. We obtained 
a continuum-subtracted Ha image, made by subtracting 
an i?-band image from a narrow-filter image at Ha (D. 
Calzetti, private communication). The pixel scale of the 
Ha image is 0.3" and the resolution is about 1". The 
Ha flux is estimated from this image by assuming the 
[Nil] /Ha ratio and fi nal Ha correction factor as given in 
iPrescott et all H20Q7D . 

As with the X-ray data, SN 2004dj was used to de- 
fine a common registration among the data sets. Un- 
saturated, foreground stars were used to register pre- 

^ http://galex.stsci.edu/GR2/ 

^ http: / / data.spitzer.caltech.edu/popular/sings/20070410_enhanced_vl /Docu 
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Fig. 7. — Identification of star-forming regions in NGC 2403. 
Top: The central 6'x6' region of NGC 2403 at FUV {GAhEX 
A1529 A). The positions of the 58 sources detected in the FUV 
band are encircled in white. The positions of the 47 sources de- 
tected at 24 /xm are encircled in black. The sizes of the circles in- 
dicate photometric apertures defined as the 3a- width of a circular 
Gaussian model fitted to the respective surface brightness distri- 
bution. The gray scale is the same as the right panel of Figure [5] 
Middle: The same region at 24 /im (Spitzer). Symbols same as for 
the FUV image. The gray scale indicate intensity levels in units 
MJy str~^ for 24 /xm image. One pixel corresponds to l"5/pixel. 
Bottom: The 19 sources common to both FUV and 24 /im images 
overlaid on the continuum-subtracted Ha image. Sixteen of the 
19 sources were also selected as Hll regions as described in § 12.41 
The gray scale indicate intensity level in the units ct s"'^ pixels ""^j 
which corresponds to 2.1 xlO~® Jy pixels"-^. The image has a 
pixel scale of 0"304. 



supernova ground-based, and GALEX images. 

3.2. Star-forming Region Definitions 

The same source finding tool applied to the X-ray im- 
ages was used to identify candidate star-forming regions 
in the FUV and 24 /im images. This tool is optimized for 
detecting point sources but can be applied to moderately- 
extended sources by increasing the characteristic size of 
the model PSF parameter. We used a circular Gaussian 
PSF in the application because most of the star-forming 
regions appeared centrally-peaked at the moderate res- 
olution of the 24 and FUV images. We examined 
their radial profiles post facto using the higher resolution 
Ha images to confirm that they are basically extended 
and either centrally-peaked or centrally-cratered so that 
the Gaussian model is adequate for estimating both the 
emission centroids and the spatial extent. 

We detected 58 bright sources in the central 6' x 6' FUV 
image (Figure [71 top panel) defined as those with a S/N 
above 10 and with a minimum of 15 source counts per 
unit uncertainty in the background. Similarly, 47 bright 
sources are detected in the 24 /im image (Figure[71 middle 
panel) with a S/N above 50 and with 100 or more elec- 
trons per unit uncertainty in the background. There are 
19 sources common to both the FUV and 24 /xm images. 
These 19 regions are shown in Figure [7] (bottom panel) 
overlaid on the Ha image. (Sixteen of these regions were 
selected as Hii regions based on their Ha brightness in 

This constitutes our star-forming region sample. The 
sample is only representative of the star-forming regions 
in NGC 2403; a larger sample could be made by includ- 
ing regions of lower S/N. Note (Figure [7]) that there 
are sources brighter than these 19 in one or the other of 
the two selection wavebands, but not in both. Examin- 
ing all available images, none of the sample regions have 
properties of foreground stars or of background AGN. 

The observed luminosities and angular sizes of the 19 
regions are listed in Table [9l We performed aperture 
photometry on the individual sources with correspond- 
ing aperture sizes for each region and waveband. Al- 
though these 19 sources are extended, in order to de- 
fine sizes of aperture, we apply the circular Gaussian 
model. The radii in each region and each band are de- 
termined as 3cr circular Gaussian widths, which contains 
99% of flux of the sources. For most regions, a sur- 
rounding annulus that extends to twice the source's aper- 
ture radius was used for background. In more crowded 
regions, a non-contiguous nearby region was used in- 
stead. Luminosities are computed using d efinitions simi- 
lar to those used by the SINGS team fe.g. JCalzetti et all 
120051 ): specifically, Ltand = ^L\ j,and for broad bands and 
Lband = SXL\j,and for uarrow bands. Luminosities were 
corrected for the Galactic e xtinction along the line of 
sig ht, Eb-v = . 040 m ag ()Schlegel et al.lll998D . with 
the lCardelh et al.l ()1989[ ) dust model. The aperture cor- 
rection values listed on the Spitzer web page were applied 
to the 24 /im data. Thumbnail images in the FUV and 
24 /im bands are shown in Figure [51 The circles shown in 
these images define the adopted radii in each band. Note 
that they are rarely concentric. Also shown in Figure [51 
are the Ha radial profiles for each region. 

The sizes of the regions in the FUV band are system- 
atically smaller than in the Ha and 24 /xm bands. The 
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Fig. 8. — 30"x30" close-up views of the 19 star-forming regions. FUV images are on the left and 24fim images are in the middle. 
Circle locations, sizes, colors, and gray scales are the same as in Figure [3 On the right are higher-resolution Hq radial profiles shown with 
best-fitting Gaussian model curves. Each radial profile extends to 5 Gaussian widths. One pixel corresponds to 0."304. (First four are 
shown here. The rest of the 15 regions are shown at the end of the paper.) 



natural interpretation is that the FUV originates from 
young star clusters whereas the source of Ha and 24 fxm 
is re-radiation from surrounding warm ionized gas and 
cold dust regions, respectively, that lie at the outskirts of 
wind-blown bubbles created by the stars. This is also the 
reason why the regions are not concentric in the differ- 
ent wavebands. The Ha radial profiles. Figure O which 
often show a crater-like morphology is also indicative of 
this shell-like structure. 



3.3. Age, Mass, & Extinction Determinations 

The five wavelength bands FUV, NUV (both from 
GALEX), 3.6 /im, 4.5 /im, and 24 /im (all from Spitzcr) 
were used to define the spectral energy distributions 
(SEDs) of light from the individual regions. These were 
compared to theoretical spectra of instantaneous star - 
bursts based o n the models of iLeitherer et al.l (|1999l ) 
as improved bv lVazguez &: Leithererj (|2005f) bv convolv- 
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TABLE 9 

Luminosities and sizes of 19 star-forming regions 





SRC 




Lfuv 
lO'^^crg s^-"- 


^FUV 
arcsec 


Lnuv 
10^^ erg s~l 


Ls.Qfim 

lO^^erg 


-^4.5^tm 

10^* erg 


L24^m 

IQSSerg 


R24^tm 

arcsec 


Lh(x 
lO^^erg 


arcsec 


10 erg s 


SF 


J073623.H 


-653649 


4.8 


± 0.3 


8.4 


2.9 ± 0.1 


4.1 


± 0.2 


2.7 ± 0.1 


4.1 ± 0.1 


8.4 


0.8 


8.5 


9.7 ± AA 


SF 


J073628.6H 


-653349 


30.2 


± 1.5 


8.3 


20.4 ± 0.6 


5.6 


± 0.2 


4.8 ± 0.1 


4.9 ± 0.1 


8.7 


7.4 


8.8 


8.9 ± 4.4 


SF 


J073635.4H 


-653507 


10.0 


± 0.5 


6.5 


5.6 ± 0.2 


5.7 


± 0.2 


3.3 ± 0.1 


2.3 ± 0.1 


7.9 


3.0 


8.0 


8.5 ± 4.2 


SF 


J073635.4H 


-653550 


4.9 


± 0.3 


7.8 


2.8 ± 0.1 


3.2 


± 0.2 


2.0 ± 0.1 


1.9 ± 0.0 


7.4 


0.7 


7.5 


1.2 ± 2.£ 


SF 


J073639.5H 


-653709 


5.4 


± 0.3 


6.7 


3.7 ± 0.1 


2.9 


it 0.2 


1.8 ± 0.1 


3.1 ± 0.1 


9.7 


1.2 


9.8 


7.7 ± 4.e 


SF 


J073641.6H 


-653806 


18.6 


± 1.0 


7.0 


10.6 ± 0.3 


6.3 


± 0.2 


4.7 ± 0.2 


5.0 ± 0.1 


9.1 


4.5 


9.2 


7.1 ± 4.2 


SF 


J073642.3H 


-653557 


4.3 


± 0.3 


6.0 


3.3 ± 0.1 


3.7 


± 0.3 


2.3 ± 0.2 


4.7 ± 0.1 


11.8 


1.0 


11.9 


21.6 ± 6.E 


SF 


J073645.4H 


-653700 


148.9 


± 7.5 


6.8 


106.5 ± 3.2 


30.8 


± 0.6 


22.7 ± 0.4 


43.9 ± 0.9 


8.5 


16.7 


8.6 


24.8 ± 6.a 


SF 


J073648.9H 


-653633 


23.1 


± 1.2 


7.3 


19.8 ± 0.6 


8.3 


± 0.4 


7.2 ± 0.2 


5.2 ± 0.1 


7.7 


1.5 


7.8 


9.0 ± 4.2 


SF 


J073652.3H 


-653646 


69.6 


± 3.5 


7.5 


53.3 ± 1.6 


19.7 


± 0.5 


13.9 ± 0.3 


36.4 ± 0.7 


9.3 


15.2 


9.4 


63.0 ± 9.2 


SF 


J073652.7H 


-653620 


7.5 


± 0.4 


6.7 


5.9 ± 0.2 


2.6 


± 0.4 


0.8 ± 0.1 


3.7 ± 0.1 


8.7 


1.1 


8.8 


7.9 ± 4.S 


SF 


J073654.H 


-653547 


13.1 


± 0.7 


8.2 


12.0 ± 0.4 


3.7 


± 0.4 


2.5 ± 0.2 


10.9 ± 0.2 


10.4 


2.0 


10.5 


24.1 ± 6.E 


SF 


J073657.8H 


-653723 


16.3 


± 0.8 


7.3 


11.2 ± 0.3 


21.3 


± 0.5 


14.5 ± 0.3 


18.2 ± 0.4 


8.2 


7.8 


8.4 


11.8 ± 5.C 


SF 


J073659.6H 


-653508 


22.6 


± 1.2 


9.0 


16.3 ± 0.5 


14.0 


± 0.5 


9.8 ± 0.3 


9.5 ± 0.2 


13.6 


4.0 


13.8 


26.7 ± 7.f 


SF 


J073700.4H 


-653524 


9.0 


± 0.5 


7.3 


7.3 ± 0.2 


2.1 


± 0.2 


2.0 ± 0.1 


4.7 ± 0.1 


11.7 


3.0 


11.9 


21.9 ± 6.E 


SF 


J073706.6H 


-653638 


150.9 


± 7.6 


9.9 


107.5 ± 3.2 


81.9 


± 1.4 


61.2 ± 1.0 


168.1 ± 3.4 


9.0 


35.7 


9.1 


98.0 ± 11.2 


SF 


J073710.4H 


-653540 


3.9 


± 0.3 


6.9 


5.2 ± 0.2 


3.5 


± 0.2 


2.5 ± 0.1 


4.3 ± 0.1 


7.8 


1.2 


7.8 


0.9 ± 2.£ 


SF 


J073717.4+653829 


18.7 


± 1.0 


6.8 


12.3 ± 0.4 


2.1 


± 0.1 


2.0 ± 0.1 


4.2 ± 0.1 


7.2 


3.2 


7.3 


0.5 ± 3.C 


SF 


J073718.0+653347 


16.9 


± 0.9 


7.4 


11.1 ± 0.3 


5.4 


± 0.2 


9.4 ± 0.2 


13.1 ± 0.3 


9.1 


6.7 


9.2 


5.2 ± AA 



Note. — FUV, NUV, and Ha luminosities are Galactic exinction corrected. 



TABLE 10 

Derived values of 19 star-forming regions 



SRC 


Mass 
10^ Mq 


Age 
Myr 




10^1 cm-i 


Lw 
lO^^erg 


r int 

lO^^erg s-1 


r int 

IQJOerg s~^ 


lO'^^erg 


SF J073623. 1+653649 


0.4 


9 


0.67 


1.7 


5.2 


9.8 ± 0.3 


1.3 


23.5 ± 10.7 


SF J073628. 6-1-653349 


0.7 


8 


0.22 


0.8 


9.9 


33.8 ± 0.3 


8.6 


13.7 ± 6.8 


SF J073635. 4-1-653507 


0.8 


16 


0.32 


1.0 


7.5 


12.8 ± 0.2 


3.7 


14.6 ± 7.2 


SF J073635. 4-1-653550 


0.3 


11 


0.44 


1.2 


3.5 


7.4 ± 0.2 


0.9 


2.2 ± 5.6 


SF J073639. 5-1-653709 


0.3 


9 


0.5 


1.4 


3.9 


8.8 ± 0.2 


1.7 


16.1 ± 9.7 


SF J073641. 6-1-653806 


0.7 


9 


0.35 


1.1 


9.0 


24.8 ± 0.3 


5.7 


12.7 ± 7.8 


SF J073642. 3-1-653557 


0.4 


9 


0.65 


1.6 


5.2 


08.5 ± 0.3 


1.6 


50.2 ± 15.1 


SF J073645. 4-1-653700 


2.4 


1 


0.3 


1.0 


26.7 


185.6 ± 0.8 


20.6 


49.2 ± 12.5 


SF J073648. 9-1-653633 


0.9 


9 


0.25 


0.9 


11.6 


26.9 ± 0.3 


1.8 


14.6 ± 6.9 


SF J073652. 3-1-653646 


1.5 


1 


0.44 


1.2 


16.7 


104.7 ± 0.6 


20.8 


142.5 ± 20.8 


SF J073652. 74-653620 


0.3 


8 


0.44 


1.2 


4.2 


11.3 ± 0.3 


1.4 


15.0 ± 8.2 


SF J073654. 1-1-653547 


0.3 


2 


0.54 


1.4 


4.7 


22.5 ± 0.4 


2.9 


61.8 ± 16.7 


SF J073657.8-(-653723 


17.5 


60 


0.67 


1.7 


0.0 


33.2 ± 0.5 


12.7 


27.5 ± 11.5 


SF J073659. 64-653508 


10 


50 


0.4 


1.2 


0.0 


32.2 ± 0.4 


5.4 


49.5 ± 14.4 


SF J073700.4-I-653524 


0.3 


7 


0.46 


1.3 


4.7 


13.8 ± 0.3 


4.2 


43.8 ± 13.0 


SF J073706. 6-1-653638 


6.4 


1 


0.71 


1.7 


71.2 


325.2 ± 1.3 


59.3 


302.5 ± 34.5 


SF J073710.4-I-653540 


0.4 


9 


0.57 


1.5 


5.2 


6.9 ± 0.3 


1.8 


2.1 ± 6.5 


SF J073717.4-I-653829 


0.2 


2 


0.27 


0.9 


3.1 


22.4 ± 0.3 


3.9 


0.8 ± 5.7 


SF J073718. 0-1-653347 


0.4 


1 


0.56 


1.5 


4.5 


29.7 ± 0.4 


9.9 


14.1 ± 11.9 



Note. — Listed A^h values include a Galactic component of A^h = 0.4 X 10^^ cm ^. 



ing these spectra with the spectral response functions 
of the appropriate filters. For computing the theoreti- 
cal spectra, we assumed the iKroupal (|2001[ ) initial mass 
function (IMF), namely a power-law IMF on the range 
0.1-100 Mq with a break at 0.5 Mq and slopes of 1.3 
below and 2.3 above the break. A solar metallicity was 
assur ned: consistent with es t imates from the lite rature 
(e.g.. iMartin k BellevI II996I : iGarnett et all 119971 ). We 
accounted approximately for the wavelength-dependent 
at tenuation by du st by using the starburst dust model 
of iCalzettil ()2001l ). We modeled the re-emission of this 
radiation in the mid-IR by assuming energy conserva- 
tion and a blackbody emission profile at a temperature 
of 75 K. This temperature was estimated by fitting a 
blackbody curve to the Spitzer MIPS 24, 70, and 160 fim 



measurements from the brightest of our star-forming re- 
gions. (This region is fully resolved and isolated even 
in the 160fim image so that a distinct local background 
could be identified.) Varying the blackbody temperature 
by ~30% has little impact on our results. However, much 
higher or lower temperatures result in unrealistic esti- 
mates of the age and mass of several of the star-forming 
regions. 

By using five-band SEDs we are thus able to sample 
the UV region which is most sensitive to age in young 
star clusters yet is most affected by dust attenuation, the 
near-IR at 3.6 and 4.5 /im that is least affected by dust 
and hence is sensitive primarily to cluster mass, and the 
mid-IR 24 /im which lies within the dust emission band 
that extends from several microns out to much longer 
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wavelengths. The 24 /im waveband is longward of the 
PAH emission features which greatly complicate the dust 
emission spectrum (e.g.. iPirez- Gonzalez et"alll2006 ') but 
shortward of the lower spatial resolution long wavelength 
Spitzer bands. 

We used a least-squares fit of the model SEDs to the 
observed luminosities using equal weights for each spec- 
tral band. Operationally, the model extinction and clus- 
ter mass parameters were allowed to vary, at a fixed clus- 
ter age, until a minimum found. We then re- 
peated this process for discrete cluster ages ranging from 
1 to 200 Myr (in steps of 1 Myr for ages <20 Myr and 
steps of 10 Myr otherwise) to obtain the overall best- 
fitting age, mass, and extinction combination. That is, 
we do not interpolate the model spectra to intermediate 
ages. 

The results are summarized in Table 1101 In addition 
to the best-fitting mass, age, and extinction values based 
on the SEDs fitting. Table [TU] also includes the corre- 
sponding hydrogen column density, Nh = (Ay/i?y)5.8 x 
10^^ cm~^, assuming an optical parameter Ry = 3.1, and 
the corresponding observed Ha, and FUV luminosities 
after correcting for extinction. 

Among the most massive young clusters in our sample 
are the gi ant Hii region s NGC 2 403-I and NG C2403-II 
studied bv lDrissen et al.l (jl999f) . IDrissen et al.l estimate 
the ages of these clusters as between 2 and 6 Myr based 
on their detection of Wolf-Rayet stars in these regions. 
A crude estimate of the initial cluster mass can be made 
from the number of WR stars detected. This estimate is 
very sensitive to the upper mass cutoff assumed of the 
initial mass function and, of course, to the true age of 
the cluster. Interestingly, the age and m ass values ob- 
tained from the studies of Drisscn et al.J are consistent 
with our finding s and indepen d ently validates our SED 
fitting method. iGarnett et al.l ()1999D estimated line-of- 
sight extinction to several of our sample star-forming 
regions from Balmer line ratios. Our estimates are in 
agreement with these results within Ay of 0.3 mag. For 
example, we derive Ay = 0.7 mag for the brightest Hii 
region, compared to Ay = 0.6 mag from Balmer line 
ratios (|Garnett et al.|[l999l ) and to Ay = 0.9 mag from 
extinction toward blue giants in the region (jPrissen et al.l 
[1991 . 

The best-fitting ages of the regions fall into three dis- 
tinct ranges, <3 Myr, from 7 to 16 Myr, and >50 Myr 
which we designate young, intermediate, and old, respec- 
tively. We believe the two regions designated as old are 
likely younger than their best-fitting ages. Inspection of 
the distribution of the fit statistic shows, for these 
two old regions, that there is a second local minimum at 
13 Myr which is probably a more realistic age estimate. 
Note that they have moderately high FUV luminosities, 
by selection, which, for their fitted age, forces their best- 
fitting mass to be very high relative to the other groups. 
There is not a similar trend in the fitted extinction val- 
ues. 

There are distinct differences among the members of 
the remaining two age groups. Members of the young 
group are typically much more luminous, per unit mass, 
at UV, Ha, and 24 /im compared to the intermediate 
age group members. Based on the star cluster model, 
the UV luminosity is expected to rise slightly during the 



first 4 Myr then drop rapidly as the most massive stars 
become supernovae. Similarly, the Ha luminosity, which 
comes from recombination of circumstellar gas photoion- 
ized by massive stars, also drops quickly as these massive 
stars disappear. The observed trends are consistent with 
this scenario. 

Dust emission, accounting for most of the 24 /im lumi- 
nosity, also represents re-radiation of starlight. It is most 
sensitive to UV but also responds to longer wavelength 
light. One would expect, therefore, that the 2Afim lu- 
minosity per unit mass would decrease more slowly than 
would the UV and Ha luminosities. However, as the 
cluster evolves and its wind-driven bubble expands, sur- 
rounding dust clouds will evaporate and the solid angle 
subtended by remaining dust clouds will decrease. This 
could account for the rapid drop with age observed in 
the 24 /im luminosity. 

These trends do not apply to the near-IR 3.6 and 
4.5 /im luminosities. Light at these wavelengths comes 
directly from stars of practically all masses (ages) and is 
not strongly affected by dust. Thus, the near-IR lumi- 
nosity per unit mass is nearly independent of the age of 
the underlying star cluster; as observed. 

One puzzling result is that the young star forming re- 
gions in our sample are much more massive than the 
intermediate-age regions. The average mass of the young 
regions is 1.9x10^ M© and only 0.5x10^ Mq for the 
intermediate-age regions. Notably, there is one region 
that is 3 times as massive as any other. Excluding this 
region, the average mass of the young group is still twice 
that of the intermediate group. This may be an indi- 
cation of a selection bias. If more massive star clus- 
ters more efficiently destroy or disperse surrounding dust 
clouds, then our selection criterion that requires regions 
be bright in both UV and 24 /im would tend to select 
against more massive clusters. Naturally, the destruc- 
tion of dust takes time so this bias may only work against 
massive intermediate-age star-forming regions. Compari- 
son of the two panels of Figure[7]shows that there are sev- 
eral regions bright in UV that lack a strong IR counter- 
part. In particular, there is a very UV-bright region just 
below the center of the field at Lpuv ~ 7 x 10**° erg s~^. 
This cluster was not selected because it lacks strong dust 
emission. It is likely of intermediate age as it also is weak 
in Ha emission. Assuming an age of 10-15 Myr suggests 
a mass in the range of about 2x10^^ to 4x10^ based on 
the values of ipuv ^'^^ the clusters listed in Table [TOl 
This makes this cluster more massive than all but the 
two most massive young clusters. 

Grude estimates of the star-formation rate in the cen- 
tral regions of NGG 2403 can be made from the val- 
ues listed in Table [TOl The total mass in stars in the 
young age group clusters is 1.1x10® M©. These stars 
were formed in the recent 2 Myr interval for a mean star- 
formation rate of 0.5 Mq yr" . Adding the contribution 
from the intermediate-age group gives 1.7x10® Mq of 
stars formed over a 16 Myr period or a 0.1 Mq yr~* av- 
erage rate. These rates can be compared to rates deduced 
from Ha luminosities. For the young age group clusters 
this rate is only 0.1 Mq yr~* and is only slightly higher 
when the intermediate-age clusters are included (see Ta.- 
ble[Tni). The two methods give quite different estimates. 
It is unlikely that the escape fraction of ionizing radia- 
tion is large enough to compensate for the discrepancy. 
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Fig. 9. — Co-added spectrum of the 6 young star-forming regions. 
The best-fitting model, an absorbed two-temperature apec model 
with the hydrogen column density allowed as a free parameter, and 
the fit residual are also shown. The spectrum has been rc-binned 
to have a width of 73 eV for display purposes. 

Finally, the values can be compared to the galaxy-wide 
star-formation rate of 1.3 Mq yr~^ based on the Ha lumi- 
no sity and assuming a nom inal A(Hq!) = 1 mag reported 
bv lKennicutt et al.l (|2003[ ). This suggests that 10% or 
up to 38% of the current galaxy-wide star formation is 
occurring within these 19 star clusters. This is not an 
unreasonable estimate since the brightest few clusters in 
the sample clearly dominate the relevant emission from 
NGC 2403. 

3.4. Comparison to X-ray Properties 

In most cases, the X-ray emission from these individ- 
ual star-forming regions is too faint to be detected as a 
discrete source using our standard source-detection al- 
gorithm. Therefore, we define the X-ray source radius 
for each region to be the largest of the values deter- 
mined from fitting the surface brightness distribution in 
the UV, mid-IR, and Ha bands. Furthermore, we com- 
bine (stack) the spectra of the young and of the inter- 
mediate age regions to accumulate sufficient counts for 
fitting of the spectra of these two age groups. For this 
purpose, we used the first observation only, in which all 
regions are imaged on the S3 CCD. We used the same 
background as was used for analysis of the residual emis- 
sion fi j2.5p and apply the same one- and two-temperature 
thermal emission models described in §2 to the stacked 
spectra. Due to low counts, only the one-temperature 
models are applied to the intermediate age region stacked 
spectrum. We calculated mass-weighted Ay from the 
SED fitting fi j3.3p for both groups and calculated a cor- 
responding A'h, namely 1.4x 10^^ cm~^ for the young and 
1.1x10^^ cm~^ for the intermediate regions. The fitting 
results are shown in Tables [TT] and [T^] for the young and 
intermediate age star-forming regions, respectively. The 
spectra with best-fitting models are shown in Figure [9] 
for the young regions and Figure [10] for the intermediate 
regions. 

Not surprisingly, the star-forming regions have X-ray 
temperatures and densities similar to the Hii regions se- 
lected by their Ha brightness (Table [6] and ^2A^ . These 
values are somewhat lower than for isolated SNRs (Ta- 
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Fig. 10. — Co-added spectrum of the 11 intermediate-age 
star-forming regions. The best-fitting model, an absorbed one- 
temperature apec model with the hydrogen column density allowed 
as a free parameter, and the fit residual are also shown. The spec- 
trum has been re-binned to have a width of 73 eV for display 
purposes. 



ble[3]) and the densities are higher than deduced for the 
residual emission (Table [5|) . The largest differences are 
found for the intermediate-age group: the derived elec- 
tron density and pressure are a factor of two lower than 
that of the Hii regions. For models allowing the interven- 
ing column density to vary, the resulting A^h are compa- 
rable to the Galactic value or about a factor of ten lower 
than for any of the other source populations considered 
including the young star-forming region group. 

These results indicate that the star clusters begin to 
show their age even after just 10—20 Myrs. Their hot 
gas densities and hence pressures begin to fall although 
their temperatures remain a moderate 2—3 MK. These 
values of density and pressure are still higher than the 
residual emission analyzed in ^2.51 Note that most of 
the LyC-producing stars will have become SNe by the 
age of 10—20 Myrs so that the traditional current star- 
formation tracer. Ha, will no longer be strong. 

We can estimate the X-ray luminosities of the individ- 
ual star-forming regions by scaling by the ratio of counts 
in the stacked spectra (for each age group) to the num- 
ber of counts detected in individual regions. The results 
are listed in Table [101 These luminosities have been cor- 
rected for extinction local to each region using the values 
of TVh tabulated in column 7 of Table [TO] (derived from 
their best-fitting Ay) and for Galactic extinction. The 
luminosities were computed using the two temperature 
model with variable Nn for the young region group, and 
the one temperature model with variable A'h for the in- 
termediate age regions. 

Figures [11] through [TO] compare these intrinsic X-ray 
luminosities to the (extinction-corrected) FUV, 24 ^m, 
and Ha luminosities, respectively. In all cases, there is 
only a weak correlation between X-ray luminosity and 
these star- formation indicators. This is due to the large 
scatter in the X-ray luminosities. 

A simple calculation shows that in one or two regions 
a relatively high X-ray luminosity may be the result of a 
faint undetected X-ray point source in the region. From 
the XLF, ^2.21 we expect about 20 X-ray point sources 
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TABLE 11 

Young regions X-ray fitting parameters 



I'it Parameters 


apec 


apec 


apec 


apec+apec 


apec+apec 


Nu (102lcm-2) 
Te, (MK) 

Normalization Ki " 
Te, (MK) 


0.4 

4+0.3 

1 9+0.2 
-^■^-0.2 


1.4 

Q+0.5 

9 q+0.4 
^■•^-0.5 


^- ' -1.1 
9 1 +0.3 
^■-^-0.5 

2g+1600 
— 28 


0.4 
9 0+0.4 
^■^-0.3 
9 q+0.4 

58+lf 


2 9+^'" 
^■^-2.0 

^•°-0.6 
^2+111 
^^-11 

6 3+1° 


Normalization K2 
erg s"! 


9 Q+0.2 
^■'-'-0.2 


2 0+°-3 
^■^-0.3 


2 0+°-^ 

Z.U_j 2 


1 1+1-^ 
9 9+0.3 
^■^-0.3 


1 0+2'5 

2 1+0.7 


L^'/IO^^ erg s-i 


2.3 


4.2 


43.7 


2.7 


12.1 


C-stat/bin 


167/108 


166/108 


154/108 


145/108 


144/108 


Derived Parameters'* 


V [/] (10''^ cm^) 


2.0 


2.0 


2.0 


2.0 


2.0 


[/"i/^] (cm-3) 


0.03 


0.04 


0.13 


0.54,0.02 


0.13,0.03 


-P/fc [/"I'^^l (10^ K cm-3) 
[Z+i/^] (105 M0) 


1.8 
0.6 


2.2 
0.8 


5.5 
3.1 


28 
1.3 


4.2 
3.5 


Eth [/+l/'] (10^3 erg) 


0.7 


0.9 


2.3 


2.7 


3.4 


tc (Myr) 


40 


21 


2.5 


90 


6.2 



<^ K = (10-y/47rD^) JneUpdV 

''volume filling factor scaling in [ ] in Column 1. 



TABLE 12 

Intermediate age regions X-ray fitting parameters 



l"it Parameters 


apec 


apec 


apec 


A^H (102icm-2) 


0.4 


1.1 


Q 4+5.1 
^•*-o.o 


Tei (MK) 


q 1+0.3 
J^-0.4 


0+0.5 
■^■U-0.5 


q 1+0.5 


Normalization Ki " 


U.D_Q 2 


^^■'-'-0.2 




Lx/lO^^ erg s"! 


^•'J-0.2 






L^'/IO^'^ erg s-i 


1.2 


1.8 


1.2 


C-stat/bin 


98.5/108 


98.7/108 


98.5/108 



Derived Parameters' 



V [/] (lO'''^' cm'^') 


4.0 


4.0 


4.0 


ne (cm-») 


0.01 


0.02 


0.01 


■P/*: [/"l^^l (lO'^ K cm-3) 


0.8 


1.0 


0.8 


(105 Mq) 


0.6 


0.8 


0.6 


Eth (1053 erg) 


0.7 


0.9 


0.7 


tc (Myr) 


60 


50 


60 



K = {10-"/4ttD'^) JuenpdV 



''volume filling factor scaling in [ ] in Column 1. 

with 0.5—2.0 keV luminosities in the range spanned by 
our sample of star-forming regions and within the central 
2'.5 radius disk of NGC 2403. Scaling this to the total 
area occupied by the star-forming regions results in an 
expected 1.7 X-ray point sources in this luminosity range 
somewhere within our sample. (Of course, the brighter 
the point source, the lower the probability of coincidence 
since the XLF scales as L~°-^^.) 

3.5. Comparison to Theoretical Expectations 

Stellar winds and supernovae within a young star clus- 
ter create a hot low-density cavity in the interstellar 
medium (ISM) which persists beyond the ~40 Myr life- 
times of the OB stars in the cluster. The cumula- 
tive effect of winds and individual SNe is to gradually 
sweep the ISM into a thin dense shell analogous to the 
shell around a stellar wind bubble. The simple physi- 
cal model of a stel lar wind bubble (jCastor et alJ 119751 : 
iWeaver et al.l Il977| ) can be taken over directly to de- 
scribe, qualitatively, these large Hii r egions or superbub- 
bles associated with OB associat ions (|McCrav fc Kafatoi 
ll987tlMac Low fc McC"ravlll988f) . In this model, massive 
stars inject kinetic energy into their surroundings creat- 



ing a freely expanding wind. Hot shocked gas surrounds 
this free wind region and occupies most of the volume of 
the bubble. Surrounding this hot gas is the dense cold 
swept-up shell. Beyond the shell is the ambient interstel- 
lar medium. Diffuse X-ray emission from star-forming 
regions arises from the tenuous shocked gas in the bub- 
ble interior and from higher density gas evaporated from 
the shell. The dense evaporated component may domi- 
nate the X-ray emission which scales as the square of the 
density. 

Although analytic versions of the original basic model 
have been applied widely in the literature, we note that 
numerous simplifying assumptions are made which may 
be poor approximations in reality. Among these are 
spherical symmetry, constant energy injection rate, adia- 
baticity, neglect of magnetic fields and of turbulence, and 
a homogeneous ambient medium. For our present pur- 
poses, the most critical of these is the assumption of a 
constant energy injection rate or mechanical luminosity, 
Lw Numerical simulations (e.g., [Strickland fc StevensI 
Il999( ) show that the expected X-ray luminosity is roughly 
proportional to and that is far from constant. 

The starburst synthesis models we have used here 
(Leitherer & Hcckman 1995) show that L„ (injected by 
winds from stars of initial mass >15 Mq and by suitably 
time-averaged SN explosions) is a function of the age of 
the star-forming region, the mass of the star cluster, and 
the assumed IMF. Specifically, only stellar winds from 
massive stars contributes to the mechanical luminosity 
during the first 3 Myrs. Then supernova explosions be- 
gin to contribute resulting in an increase in by about 
a factor of 4 from an age of 3 to 6 Myrs. The lumi- 
nosity decreases again as the contribution from stellar 
winds declines due to the decreasing number of O stars. 
From '^10 Myrs and up to the last supernova, which is 
around 40 Myrs, is almost constant then finally falls 
quickly to very low values. The young and intermedi- 
ate star-forming regions studied here have ages in the 
range where changes significantly on short timescales 
so that our estimates of L„ (Table [T0|) are reliable only 
to about a factor of four. 

Figure [13] displays the (extinction-corrected) X-ray lu- 
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Fig. 11. — Lpuy vs. Lx for the 19 star-forming regions. Lumi- 
nosities are corrected for intrinsic and Galactic extinction. Trian- 
gles represent young (less than 2 Myrs), circles intermediate (2—20 
Myes), and squares old (above 20 Myrs) star- forming regions. The 
straight line indicates the best-fitting linear function with a slope 
of (5.5±0.8)xl0-5. 




10 100 



Fig. 12. — L24(im vs. Lx for the 19 star-forming regions. Lumi- 
nosities are corrected for intri nsic and Galactic extinction. Sym- 
bols are the same as in Figure [TT] The straight line indicates the 
best-fitting linear function with a slope of (1.9±0.2) X 10"*. 
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Fig. 13. — vs. Lx for the 19 star-forming regions. Luminosi- 
ties are corrected for intr insic and Galactic extinction. Symbols are 
the same as in Figure lTTI The straight line indicates the best-fitting 
linear function with a slope of (3.5±0.5) X 10~^. 



minosity estimated in N3.4I against the mechanical lumi- 
nosity computed from the (jLeitherer fc Heckmanlll995f) 
models using the age and mass estimates of the star- 
forming regions deduced in §3.31 (Table [T0|) . Filled tri- 
angles depict young regions, and filled circles represent 
intermediate-age regions. We fitted with linear func- 
tion for each group, and obtained the slope Lx/L^. 
The resulting slopes are 0.0033±0.0005 for the young 
regions, and 0.0017±0.0005 for the intermediate-age re- 
gions. This is consistent with the Chandra result of 
the second brightest Hii region in the lo cal group, NGC 
604, w hose Lx/L^j is estimated as 0.002 (jTiillmann et al.l 

\2om . 

Numerical simulations (jStrickland fc Stevens! Il999f) 
found that Lx/L^ for a thin galactic disk model is 
0.5—2%, using an Lx appropriate for the ROSAT band. 
This must be adjusted for the different sensitivity of 
Chandra compared to ROSAT for the observed spectral 
shape. Here, the characteristic electron temperature is 
2—3 MK so that the fraction of the luminosity for Chan- 
dra band (0.5-2.0 keV) to ROSAT band (0.1-2.5 keV) 
is 0.5 to 0.6. This translates Lx/L^ to 0.25 to 1% for 
the Chandra band which is consistent with our results. 

Figure [15] displays the (extinction-corrected) X-ray 
emission against the estimated cluster age. Here, the 
X-ray luminosities are scaled to a fixed cluster mass 
of 10^ M0 for easy comparison to Lw(0- The solid 
curve represents the best-fitting scaling, Lx = O.OOlTiw, 
from Figure [Til for the intermediate-age clusters and the 
dashed line is the factor-of-two higher scaling appropriate 
for the younger clusters. Note the roughly factor-of-four 
change in L^ at i ~ 4 Myr as described above. Unfortu- 
nately, there is a considerable amount of scatter among 
the individual regions. For this reason, the LxjL^ ratio 
is not a sensitive independent measure of the age of a 
cluster. 

4. DISCUSSION OF INDIVIDUAL STAR-FORMING 
REGIONS 

We have examined 19 star-forming regions in the cen- 
tral part of NGC 2403. Their basic properties are listed 
in Tables [51 and [TPl They range in age from 1 to 16 Myr 
with two much older regions at 50 and 60 Myr. These 
latter ages are likely poor estimates as these two regions 
are also strong Ha sources which argues against an old 
age. The remaining star- forming regions fall into distinct 
age groups; young (1-3 Myr) and intermediate (7 to 16 
Myr). The bulk X-ray properties of the hot gas within 
these two age groups, from model fits to their stacked 
X-ray spectra, are summarized in Tables [TT] and 1121 re- 
spectively. 

There is very little change in the X-ray temperatures 
and very little change in the X-ray luminosity per unit 
star cluster mass with age (Figure [T5[). There is also no 
significant evolution in the ratio of the X-ray luminos- 
ity to the mechanical luminosity (Figure [T^ . Thus, the 
basic X-ray observables, temperature and luminosity, do 
not trace evolutionary changes for this sample of young 
star- forming regions. 

However, the density of the hot gas derived from the 
spectral fitting and hence the gas pressure in the young 
star clusters is about a factor of 2 higher than in the inter- 
mediate age clusters. This is consistent with the trends 
deduced previously (^2]) where the youngest sources, the 
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Fig. 14. — Mechanical luminosity, Lw calculated from Star- 
burst99 scaled by the mass and age of each region is shown against 
the int rinsic X-ray luminosity. Symbols are the same as in Fig- 
ure 111! (the two regions older than 20 Myr are omitted because 
they have very low mechanical luminosities according to Star- 
burst99). The Unes indicate best-fitting slopes of (3.4±0.5) X lO""^ 
and (1.7±0.5) X 10^'' fitted to the young and intermediate-age re- 
gions, respectively. 
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Age 

Fig. 15. — Intrinsic X-ray luminosity is shown against age of 
star forming regions. Luminosites are scaled to an initial mass of 
10^ M0. Symbols are the same as in Figure [TT] Open diamonds 
depict the average for young and intermediate regions. The solid 
curve indicates the X-ray luminosity scaled from the mechanical 
luminosity, L„, assuming Lx/iw = 0.0017. The dashed curve 
indicates the X-ray luminosity scaled from L„ assuming Lx/iw = 
0.0034. 

SNRs, have the highest densities and pressures foUowed 
by the (Ha-selected) Hii regions with the residual X-ray 
emission having the lowest density and pressure. 

There are three plausible scenarios for the nature of the 
hot residual gas consistent with these observed trends. 

(1) If the decrease in density is temporal, then the resid- 
ual emission comes from gas initially heated at times 
considerably more than ^^20 Myr in the past (the age of 
the intermediate group of star-forming regions); i.e., it is 
a relic of past star-formation activity. This is consistent 
with the absence of a spatial correlation between trac- 
ers of recent star formation and this residual emission. 

(2) The hot residual gas may be escaping from active 
star-forming regions into lower-density voids in the disk 



ISM. These low-density voids may be, for instance, lo- 
calized remnants of past star-forming activity. If this 
were the case, then we would expect some of this hot 
gas to be surrounding the star-forming regions analyzed 
here. We checked this possibility by taking successively 
larger X-ray source sizes for the 19 star-forming regions 
to estimate the true source extent. We found that for 
sizes up to twice the values adopted above (which were 
the maximum radii estimated from the UV and mid-IR 
images) there is a clear net X-ray excess (above the back- 
ground). Nevertheless, on larger scales, the morphology 
of the residual hot gas does not correlate well with the 
star-formation tracers. (3) The hot residual gas may have 
moved out of the plane of the galaxy down the density 
gradient into the halo. However, again, as there is no 
evidence from the distribution of the residual gas that 
strongly correlates it with star- forming regions, it must 
be a relic of past activity. In fact, there are no regions of 
active star formation that have the signature of blowout 
from the disk. Blowout requires that the bubble size be 
roughly a few times the density scale height. The largest 
star-forming regions in the central regions of NGC 2403 
are ~200 pc compared to a canonical scale height for 
la te-type spirals of 250—500 pc (see also the discussion 
in [Strickland et all 2004). 

We note that the age of the hot residual gas 
inferred from all of these scenarios is roughly the 
same as the timescale for gravitational instabilities to 
cause shells sur rounding OB associations to fragment 
(jMcCrav fc Kafatos T987i) at which point hot gas can 
vent out into the surroundings. It is also an age at which 
strong ionizing continua from massive stars has ended 
though SNe may still be active. The current cooling 
time for the residual gas exceeds 20 Myr and may be 
as high as 340 Myrs. Thus, regardless of the actual lo- 
cation of the residual gas (disk or halo) and the time 
elapsed since it formed (<10— 20 Myr if escaping from 
active star-forming regions or perhaps factors of a few 
longer if remnants of past activity), this gas is likely to 
continue to radiate at low levels for a time longer than 
the characteristic timescales for localized star formation. 

5. SUMMARY 

We have revisited the study of the current and recent 
star formation activity in the central regions of normal 
late-type spiral galaxy NGC 2403 by re-evaluating the ex- 
isting X-ray data and including supporting observations 
at other key wavebands. These include Spitzer images 
that trace dust heated by massive stars, GALEX im- 
ages showing the young stars, and Ha observations that 
reveal nebulae photoionized by massive stars and SNe. 
Through analysis of this multi-wavelength data, we have 
obtained estimates of the mass, age, and line-of-sight ex- 
tinction towards numerous young star clusters. With few 
exceptions, these clusters and their environs are weak X- 
ray sources. Even the most powerful, a giant Hii region 
comparable to 30 Dor in the LMC and NGC 604 in M33, 
radiates at only ^^.3x10'^^ erg s~^ in the Chandra X-ray 
band which is about 0.5% of its estimated mechanical 
luminosity. 

We have also shown that, after carefully accounting for 
the point-source (X-ray binary) population, SNRs and 
bright Hii regions, there remains a residual X-ray com- 
ponent pervading the central few kpc of NGC 2403. This 
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component cannot be accounted for by faint sources be- 
low our detection limit. It is likely diffuse hot gas but not 
strongly correlated with current star-formation activity. 
It is likely a relic of star formation activity occurring 
some 20 Myrs (the age of the intermediate age regions) 
or older. 

The geometry of the residual hot gas is not well de- 
termined, which leads to additional uncertainties in its 
physical properties. It is not clear whether this gas is 
confined to the disk of the galaxy or resides in the halo. 
Nonetheless, it seems to be more centrally-located than 
the star-forming regions. Since the amount of hot gas 
represented by this residual emission is substantial, it im- 
plies that star-formation activity was much higher in the 
past. Combined with its central loc ation, this conclusion 
is con sistent with the suggestion bv lDavidge fc Courteaul 
([2002) who suggested that an earlier episode of star for- 
mation occurred in the central region of NGC 2403 that 
has now propagated outward in the disk to its present 
radius of about a kpc from the galactic center. 
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Fig. 8 CONT.. — 30"x30" close-up views of the rest 15 of 19 star-forming regions. FUV images are on the left and 24 /xm images are in 
the middle. Circle locations, sizes, colors, and gray scales are the same as in Figure[7] On the right are higher-resolution Ha radial profiles 
shown with best-fitting Gaussian model curves. Each radial profile extends to 5 Gaussian widths. One pixel corresponds to 0."304. 



